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Abstract 

We compute the master integrals containing one massive propagator en- 
tering the two-loop electroweak form factor, i.e. the process // X, where 
// is an on-shell massless fermion pair and X is a singlet particle under 
SU{2)l X U{1)y, such as a virtual gluon or an hypothetical Z'. The method 
used is that of the differential equation in the evolution variable x = —sjm^, 
where s is the cm. energy squared and m is the mass of the W or Z bosons 
(assumed to be degenerate). The 1/e poles and the finite parts are com- 
puted exactly in terms of one-dimensional harmonic polylogarithms of the 
variable x, H{w]x), with e = 2 — D/2 and D the space-time dimension. We 
present large-momentum expansions of the master integrals, i.e. expansions 
for |s| 3> m^, which are relevant for the study of infrared properties of the 
Standard Model. We also derive small-momentum expansions of the master 
integrals, i.e. expansions in the region |s| <C m^, related to the threshold 
behaviour of the form factor (soft probe). Comparison with previous results 
in the literature is performed finding complete agreement. 
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1 Introduction 



Present day and planned high-energy experiments often require accurate perturba- 
tive calculations. While in quantum-choromodynamics (QCD) higher order compu- 
tations are available for various processes, in the electroweak (EW) case progress is 
more difficult because of the presence of the large, non negligible, masses of W, Z 
and Higgs bosons. A full calculation of two-loop electroweak corrections to a typical 
process such as e~^e~ jX^ jj," is still outside present-day technology. It involves 
indeed the computation of a large number of massive two-loop box diagrams having 
thresholds at s = 0, m?, 4m^ and s = 9m^, assuming m ~ niw ~ mz ~ rriH- To 
have an explicit estimate of the size of two-loop electroweak corrections, we consider 
then the model process: 

/(Pi) + /(P2)-X(g), (1) 

where / (/) is a massless fermion (antifermion) on its mass-shell, Pi = = 0), 
and X is any particle without electromagnetic or weak charge, i.e. a singlet under 
the electroweak gauge group SU{2)l x U{1)y, such as a (time-like) gluon or an 
hypothetical Z'. Eq. (P) also represents the simplest process containing electroweak 
double logarithms, considered for the first time, as far as we know, in Reaction 

is the electroweak analog of the QCD processes involving a color-neutral probe, 
such as Drell-Yan production of vector bosons or muon pairs: qq^'^ W, Z, fi^fi^ 
or Higgs production by gluon-gluon fusion: gg — > H. 

At two-loop level, the annihilation in Eq. (Q) involves the emission of a pair of 
virtual bosons among 7, W, Z and if's. In the case of double photon emission, no 
mass is present in the loops and we recover the known results for the QED or QCD 
amplitudes. The latter, for dimensional reasons, are of the form: 

(massless amp.) = G{e) i-^j s^""'' = a^""-* (^j G(e)x^-"'*-2^ (2) 

where s is the cm. energy squared, G{e) is a function having a Laurent expansion in 
e = 2 — D/2 (with D the space-time dimension), n is the mass-scale of Dimensional 
Regularization and is the number of scalar propagators. We have defined: 

X = — ^ , (3) 

where m is the mass of the W, Z or Higgs bosons'*. On the other hand, if W, Z or 
i7's are also emitted, the dependence on x is not factorized in a single power and 
the amplitude has the more general form 

(massive amp.) = a* "'^ I — I G{x;e). (4) 

V a / 

The amplitudes involving a single massive propagator, i.e. a single W or Z emis- 
sion, have thresholds only at s = and s = m^. As a consequence of this simple 



The minus sign in the definition of x is inserted for convenience. 
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structure, these amplitudes can be exactly computed in terms of one-dimensional 
harmonic polylogarithms |21 Ej (a generalization of the Nielsen's polylogarithms 
lU E]) of argument x, H{w,x), with weight up to w = 4 included. In this paper 
we present the evaluation of the master integrals entering the above-mentioned am- 
plitudes, leaving to a future work the more complicated amplitudes with two and 
three massive propagators, involving also the threshold s = Am?. Since external 
fermions are massless and are taken on their mass-shell, infrared divergences (soft 
and coUinear) occur in the amplitudes. When an intermediate vector boson is ra- 
diated, its mass m acts as an infrared regulator and mass logarithms do appear in 
four dimensions, i.e. terms of the form log'^(— s/m^) with k up to two included. 
Infrared divergences related to photon emission instead are not regulated by any 
physical parameter. Our computation is done within the Dimensional Regulariza- 
tion scheme jH], in which the latter manifest themselves as poles in 1/e. These poles 
are unphysical and, in the physical cross-section, are canceled by the corresponding 
ones appearing in the real photon emission contributions or are factorized into QED 
structure functions. Ultraviolet divergences — related to short-distance behaviour 
and leading to additional 1/e poles — are also regulated within Dimensional Regu- 
larization and can be subtracted with ordinary renormalization prescriptions, such 
as, for example, the MS scheme. 

The computation of the form factor is naturally divided in two steps. 

The first one involves the reduction of the original tensor diagrams, given in 
Figs.nand|2l into a minimal set of scalar amplitudes, the so-called master integrals. 
This step is discussed in section |^ and involves in turn three different operations: 
a) the decomposition of tensor amplitudes into scalar amplitudes by means of pro- 
jectors on invariant form factors, see section 2.1; h) the transformation to linearly 
independent scalar amplitudes by means of scalar product rotations, as discussed in 
section 2.2; c) the reduction of all the independent scalar amplitudes generated to a 
small subset of them — the so called master integrals — by means of the integration- 
by-parts identities [7j and other symmetry relations, as discussed in section 2.3. 

The second step involves the explicit evaluation of the master integrals, which 
we perform in all the massive cases with the method of differential equations on the 
external kinematical invariants jH[ El 1^ 1^ 1^ ■ A set of special functions of 
the polilogarithm kind is needed to represent the pole and finite parts of the master 
integrals. We choose the basis of the (one-dimensional) harmonic polylogarithms 
13 EI • This step is described in section El 

In section|3]we list the results for the pole and finite parts of the master integrals 
containing up to five denominators included. In section El we give the expressions 
of all the scalar integrals containing six denominators. The planar topologies are 
reducible while the crossed ladders — two amplitudes — are master integrals. We 
decided to list them together because all the six-denominator amplitudes are inter- 
esting by themselves for reference. We also compare with previous results in the 
literature. Section El is devoted to the large momentum expansion |s| ^ m? of all 
the six-denominator amplitudes. As already said, this expansion is relevant for the 
study of asymptotic properties of multiple electroweak radiation. Section [71 deals 
with the small momentum expansion — i.e. the expansion for |s| ^ iv? — of all the 
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six-denominator amplitudes. In the conclusions we summarize the results obtained 
and we discuss future perspectives. To not bore the reader with too many technical 
details or formulas — which are anyway worth reporting — we added a few appen- 
dices. Appendix ^ contains the expressions of all the denominators entering the 
two- loop amplitudes. In appendix |ni we give the results of the one-loop amplitudes 
which enter the factorized two-loop diagrams, i.e. the diagrams in which the two 
loops do not have common propagators. Finally, in appendix O we list the results 
for all the independent topologies of Figs. El and [7| which are not master integrals 
but can be useful for general reference. 

2 The reduction to master integrals 

All the two-loop vertex diagrams containing six denominators with at most one 
massive propagator are reported in Figs. Hand|2P. The graphical conventions we 
use are the following. A massless propagator is represented by a wavy internal 
line, while a massive propagator is denoted by a straight segment. An external 
line carrying the light-cone momentum pi or p2 is represented by an external wavy 
line, while the external particle X carrying the general momentum q is depicted as 
a straight external line. This set of conventions is clearly a generalization of the 
well-known QED ones. In Fig. Q we consider the basic three topologies with the 
various mass configurations: 

• (a), (b), (c) and (d) are ladder diagrams, i.e. corrections to the hard vertex of 
the basic one-loop amplitude; 

• (e) and (f) are crossed ladders — the non-planar topology; 

• (g), (h) and (i) finally are the vertex-insertion diagrams, containing a correc- 
tion to the vertices with momentum pi or p2. 

The number of internal lines with different momenta entering a diagram is an indica- 
tor of the topology of the diagram. The diagrams in Fig. ^ are real six-denominator 
amplitudes and the related topologies are given in Fig. HI 

In Fig. 121 the diagrams containing self-energy corrections to the internal lines of 
the one-loop vertices are plotted. Note that, in contrast to Fig. [TJ diagrams (a), (b), 
(c) and (d) in Fig. |21 involve only five distinct denominators — one of them being 
squared — as two internal lines carry the same momentum. These diagrams then, 
do not represent real six-denominator topologies and are more naturally grouped 
among the five-denominator ones. The topologies of diagrams (a), (b), (c) and (d) 
are shown in Fig. in (h), (e), (f) and (s) respectively. Diagram (e) of Fig. |21 
involves two internal lines with the same momenta but with different masses. It can 
be expressed as a superposition of integrals belonging to topologies (s) and (t) of 

^Actually, diagram (/) in Fig. [^contains only five denominators, but we list it in this figure as 
it belongs to the same class of the six-denominator ones. 
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Figure 1: two- loop vertex diagrams with up to one massive propagator. They 
involve the ladder, crossed-ladder and vertex-insertion topologies and are real six- 
denominator topologies (see text). 
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Figure 2: two-loop vertex diagrams with up to one massive propagator involving a 
self-energy correction to the basic one-loop vertices. The related topologies are five 
and four denominator topologies (see text). 
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Fig. 13 by means of partial fractioning in the loop momentum ki. 

1 _ 1 
kf {kf + a) a 

The scalar product of two vectors is defined as: a ■ b = — ao&o + a ■ b ^. In graphical 
form, the decomposition reads: 

Finally, diagram (f) is actually a four-denominator topology, indicated in Fig. EKg)- 
To conclude, let us note that the topologies related to diagrams of Fig. |21 can 
be obtained by properly shrink one internal line (or two internal lines in the case of 
diagram (f)) among those in the diagrams of Fig. EJ 



(5) 



(6) 



2.1 Decomposition of a tensor amplitude into scalar ampli- 
tudes 

In this section we consider the decomposition of a general tensor amplitude among 
those shown in Figs. nand|21into scalar amplitudes. Even though the decomposition 
is completely general, let us pick up, for clarity's sake, a specific case: the ladder 
diagram with a mass on the external boson line (the diagram (b) in Fig. P). Omitting 
an overall constant dependent on the couplings of the theory, this diagram involves a 
double integral over the loop momenta ki and k2 and can be written, in all generality 

7. 



as 



J -f^l-r^2-f3-' 4-' S-T 6 

where {d^k} is the loop integration measure (see section H}. According to our 
routing of the loop momenta, we have defined the denominators: 

Pi = kl + a, (8) 
P2 = {Pi-k{)\ (9) 

Ps = {P2 + h)\ (10) 

Pa = A;2^ (11) 
A = {Pi~h + k2)\ (12) 
= {p2 + ki-k2f. (13) 

The numerator M^vij... depends on the spin of the particles entering the diagram and 
on the type of interaction; the indices ixvij... denote collectively Lorentz indices, 

^Note a difference of sign with respect to the more common definition. 

''^Depending on the gauge choice for vector particles, some of the denominators Pi may actually 
appear to power two instead of one: the reduction scheme is clearly unchanged. We do not consider 
the non covariant gauges, leading to additional denominators and scalar products involving a 
constant gauge vector c'^. 
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Dirac indices, etc. By using projectors on invariant form factors, we can always 
restrict ourselves to consider scalar numerators (see for example ^|). The latter 
typically involve a trace over spinor indices and a sum over the polarizations of the 
vector particles. Therefore, from now on, let us consider: 



Af{Si, S2, 5*3, 5*4, 5*5, 5*6, 5*7 
-Pl -P2 -P3 -P4 -P5 -Pe 



(14) 



where we have explicitly indicated the dependence on the seven invariants, depen- 
dent on the loop momenta ki and k2'- 

(15) 
(16) 
(17) 
(18) 
(19) 
(20) 
(21) 

The scalar numerator can be expanded in a power series of the invariants as: 

Af{Si, S2, S3, S4, 5*5, 5*6, 5*7) = b{li, l2, - ■ ■ , h) S[^ S2 S^3^ S'^ S^^^ Sq^ S!j^ , (22) 

li,--h>0 

where the (known) coefficients b{li, h-, - ■ ■ ih) also depend on the external invariant 
s, the mass squared a = m? and the space-time dimension D. With the typical 
gauge interactions, we expect the indices U to be restricted to U < 3. 



Si 


- 

— AL]^, 




S2 


- k^ 




S3 


= fci- 


Pu 


s. 


= h- 




S5 


= k2- 


Pu 


Se 


= k2- 


P2, 


Sr 


= h- 


k2. 



2.2 Decomposition into independent scalar amplitudes 

The scalar amplitudes in Eq. (fT^ . with J\f given by the expansion in Eq. are 
not, in general, linearly independent on each other; we can indeed express six of the 
invariants listed in Eqs. (jl5fl^T|) in terms of the denominators Pj. The next step is 
the reduction to a set of linearly independent scalar amplitudes in order to progress 
in the evaluation of the tensor amplitude of Eq. (|7j). 

In this work, we used two different reduction schemes of the scalar amplitudes, 
using consistency of the results as a crossed check. The first one is the scheme of 
the auxiliary denominator (or auxiliary diagram) , which is essentially based on the 
reduction of all the amplitudes to scalar ones containing formally only denominators; 
the second one is the scheme of the shifts, in which the independent scalar amplitudes 
are integrals with a subset of the scalar products at the numerator. 

In the following two sections we will recall briefiy both schemes, paying more 
attention on the first one, not discussed often in the literature, and just sketching 
the second one, for which we refer to the bibliography. Our results are presented in 
section m in the framework of this second scheme. 
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2.2.1 Auxiliary-denominator scheme 



As already said in the previous section, we aim at expressing the invariants listed in 
Eqs. ()15fl^T|) in terms of the denominators Pi appearing in the diagram. Since there 
are seven invariants and only six denominators, we introduce a seventh auxiliary 
denominator, independent from the others, in such a way to form a basis. The 
choice is of course, not unique; let us take for example: 



Pi = ih - k^f 



We then have the linear relations: 



7 

Si = ^ ^ -^ij (Pj '^1 



,7) 



(23) 



(24) 



where aj = Sj^ia and Aij an invertible 7x7 constant matrix. The non-zero matrix 
elements are given by: 



Si 

S3 

s, 

S5 

s, 

S7 



Pi -a, 

A, 
1 



-P2 
-P7 



Pi -a), 
Pi + a), 
P5-P2 



Pi -a) 
+Pr -Pe + P3-Pi + a) 
-Pj + P4 + Pi - a) . 



(25) 
(26) 

(27) 
(28) 
(29) 
(30) 
(31) 



Let us note that one can construct a fictitious Feynman diagram having all the seven 
denominators, called the auxiliary diagram |16|: 



(aux. diag.) 



Pi P2 P3 Pi P5 P& Pj 



(32) 



In our example, the auxiliary diagram is a planar double box with final momenta 
equal to the initial ones pi and i-e. a four-point forward amplitude: see Fig. 01 
In general, an auxiliary diagram has a number of independent denominators equal 
to the number of the invariants (seven in our case). The scalar numerator M can 
now be expanded in powers of the basis "vectors" {Pi}i=ij as: 



c{si,S2,---,sr)Pl^P^^---Pr 



Sl,---S7>0 



with the new coefficients c(si, S2, ■ ■ ■ , S7) being linear combinations of the old ones 
b{li, h - ■ ■ Inserting the above expansion in the expression (fTH) for the diagram. 
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Figure 3: Double box auxiliary diagram for the case of the vertex diagram in Fig-Q] 
(b). The dashed line represents the auxiliary denominator Pj. 



we obtain: 

^ \- r ^ f {d^hUd^h} 

si,---S7>0 1 2 3 4 5 6 7 

Eq. ()33|) is the main result of this section. It says that the general scalar diagram 
X is a linear combination — with known coefficients — of independent scalar am- 
plitudes involving the original denominators {Pj}j=i,6 and an additional, fictitious, 
denominator P7. A compact notation for the auxiliary scalar amplitudes is the 
following: 



Topo(ni,n2,n3,n4,n5,n6;n7) = y 
according to which Eq. is written as: 



{d^ki}{d^k2} 



pn2 pna p"4 pris pm pm ' 
-'2 -'3 -'4 -'5 -'6 -'7 



(34) 



^iPl,P2)= 5Z C?('^l, "■2, ••■'^7)Topo(ni, 712, 713,714, 715,716; 77,7). (35) 
ni , --ng<l, n7<0 

A few remarks are in order. First, the auxiliary denominator P7, if it appears, it 
only does at the numerator: it cannot clearly be generated at the denominator by 
tensor decomposition. In other words, rii = 1, 0, —1, —2, —3, ... for i = 1, 6 while 
717 = 0, — 1, — 2, — 3, ... The second remark concerns the so-called topologies. As 
already anticipated, an important quantity for a scalar amplitude is the number of 
distinct denominators with positive indices, rij > 0, which we call t: 



5^0(71,), (36) 



i=l 



with 6 {x) the unit step function defined in such a way that ^ (0) = 0. The scalar 
amplitude with ni = 2, 712 = ■ ■ ■ = uq = 1 and nj = —1 for example has 
t = 6, the maximal value for two-loop vertex functions. If a scalar amplitude 
Topo (711, 712, 713, 714, 715, 7i6; 717) has 71, < 0, the i*'^ denominator eventually appears 
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only at the numerator and the i^^ internal line is shrunk to a point: the diagram 
then has less than six internal lines. Here we give an example of a decomposition 
into auxiliary amplitudes with t = 6 and t = 5: 

y"{d%}{d%}-p^^^g^ = -1topo(1, 1,1,1,1,1;-!) 

-|Topo(l,l,l,l,l,l;0) 
+^Topo(0,l,l,l,l,l;0) 
+iTopo(l,l,l,0,l,l;0). (37) 

In general, a diagram X is decomposed into auxiliary amplitudes with decreasing 
t starting from t = 6 included, i.e. t = 6,5, 4, 3*. 

2.2.2 Shift scheme 

We described above the scheme of the auxiliary diagram to reduce a general scalar 
amplitude. Let us now present another scheme of reduction which has also been 
used in the past (see [T71 IT^ E])- The main results of this work are presented in 
sections m and El according to this second scheme, so let us briefly overview it. 

In essence, it uses shifts on the loop momenta in order to maximally simplify the 
structure of the sub-topologies. Let us consider a routing with two internal lines 
having the loop momenta ki and k2, as for instance Pi and P4 in the discussion 
above. We choose for example to keep the invariant 6*7 and we simplify, i.e. cancel, 
the factors Si{i = 1, ... ,6) appearing at the numerator with the denominators by 
means of the formulas: 



Si 


= Pi - a 


(38) 


S2 


= P4 


(39) 




= l(-P2 + Pi + a) 


(40) 




= l(+P3-P, + a) 


(41) 




= ^(+P,-P,-P,) + Sj 


(42) 




= ^{-P, + P, + Ps)-St. 


(43) 



The simplification of a scalar product with the selected denominator brings to the 
appearance of sub-topologies in the case the denominator appears to power one. In 
principle, a diagram of topology number t can have t\ different sub-topologies. In 
our case, however, some sub-topologies actually vanish because all the propagators, 
except possibly one, are massless. Furthermore, eventual discrete symmetries of the 

^Amplitudes with t < 3 vanish, in our case, as they contain a massless tadpole. 
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diagram, such as an up-down symmetry, further reduce the number of independent 
sub-topologies. 

After the cancellations between the Si and Pj have been done, we end up with 
a set of subdiagrams in which the denominators Pi and P4 containing the loop 
momenta ki and k2 may be absent. In this case, we perform shifts on the loop 
momenta to reproduce denominators containing again ki and /c2- For example, if Pi 
is absent but P2 is still present in the given amplitude, we may do the shift 

ki^pi-ki. (44) 

If Pi and P2 are absent but P3 is still present, we may do the shift ki — > ki — p2, 
and so on. In general, we construct a table of shifts on the loop momenta. With 
this step, new denominators are generated. Under (jl^ for example, 

P3 {Pi+P2~kif, 
P5 ^ (k^ + k^y, 

^ iPi+P2-h-k2Y, (45) 

denominators which were not initially present. We then extend the set of cancellation 
rules including those between the new denominators and the scalar products. The 
scalar product ki ■ ^2 cancels against the denominator [ki + ^2)^ — no arbitrareness 
— while {pi+p2 — ki)^ may he cancelled for example withp2-^i and (pi+P2 — ^i — ^2)^ 
with P2 ■ k2. The remaining scalar products are then pi ■ ki and pi ■ /c2- Depending 
on the topology number t of the diagram under consideration, we can re-express 
t of the 7 scalar products in terms of denominators and simplify such expressions 
with a power of the corresponding denominator. It follows that a diagram with 
topology number t can have at most t = 7 — t independent scalar products in the 
numerator. A six-denominator amplitude {t = 6), for example, has at most one 
scalar product, while a sunrise diagram {t = 3) has at most four distinct scalar 
products in the numerator. After all the second-step cancellations have been done, 
we go for a second step of shifts of the loop momenta, and so on. We repeat 
iteratively the cancellation-shift procedure until no more simplifications are feasible. 
The whole set of required denominators is listed in appendix^ In our computation 
we found 9 independent scalar amplitudes with six denominators, i.e. t = 6, (Fig. 
0}, 22 independent amplitudes with five denominators, (Fig. Ej), 17 independent 
amplitudes with four denominators (Fig. ^ and 4 independent amplitudes with 3 
denominators (Fig. [7j). This makes a total of 44 independent scalar amplitudes. 

Within this method, for a given topology, i.e. for a given set of denominators 
Vi^Vi,^...Vi^ (see appendix the independent scalar integrals are of the form: 

/5'*i . . . 5'** 
{d'^ki}{d^k2} _ _ X ' (46) 
11 it 

where t = 7 — t denotes the number of independent scalar products 5*4^ ■ ■ ■ Si-. Note 
that Si, ■ ■ ■ , Sj > while ni, ■ ■ ■ , nj > 1. 
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(a) (b) (c) (d) (e) 




(f) (g) (h) (i) 



Figure 4: The set of 9 independent 6-denominator topologies. 



The advantages of the auxihary diagram method are the simphcity and a shorter 
code for its implementation. The main advantage of the shift method is that a given 
topology appears in a unique form, so the CPU time needed is less. 

Let us conclude this section with a general remark. A straightforward strategy 
to evaluate X is that of computing all the independent scalar amplitudes entering 
its decomposition. In simple calculations, this is probably the simplest thing to 
do. In more complicated cases, however, there may be a large number of such 
scalar amplitudes and the individual computation of all of them may become rather 
laborious. In these cases, it is more convenient to use a set of relations connecting 
different scalar amplitudes, and to compute directly only a subset of them. The 
latter is our strategy that is described in the next section. 

2.3 Relation between independent scalar amplitudes: the 
integration by parts identities 

As we have shown in the previous section, the task of evaluating a two-loop vertex 
diagram with an arbitrary scalar numerator M is shifted to that of computing a 
set of independent scalar amplitudes either containing an auxiliary denominator: 
Topo(ni, n2, n3, n4, ns, ng; ny), or containing a selected basis of denominators and 
scalar products: Topo(ni, ■ ■ ■ , n^; Si, ■ ■ ■ , sj). Let us discuss the simpler case of the 
auxiliary diagram first. The allowed range of the indices are: ni, ■ ■ ■ ,n6 < 1 and 
rij < Q (see Eq. ((Ml)), but in practice, we expect the relevant nj's to be restricted 
to rii > —3, which however gives still a rather large set of scalar amplitudes to 
compute. Relations among different scalar amplitudes are obtained by considering 
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(a) (b) (c) (d) (e) 

(f) (g) (h) (i) (j) 




(k) (1) (m) (n) (o) 




(p) (q) 



Figure 6: The set of 17 independent 4-denominator topologies. 




(a) (b) (c) (d) 

Figure 7: The set of 4 independent 3-denominator topologies. 
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the so-called integration- by-parts (ibp) identities [7j: 



^ — J {d ki}{d ^ pni pn2 pn3 pn4 pns pne pnr^ y (47) 

with i = 1,2 and v = ki, k2,Pi,P2- The integral above vanishes because of the 
following argument. It is the space integral in D dimensions of the total divergence 
of the vector in curly brackets. It can then be transformed into the flux integral of 
this vector over a spherical surface with infinite radius, r = oo. For small enough 
dimension D, the surface integral vanishes for r — >■ cxd. It can then be defined to 
identically vanish for any D. By explicitly performing the derivations in Eq. (jTTj) 
and expressing the scalar products in terms of the assumed basis {Pi}i=ij according 
to Eqs. (j24j) ■ we obtain identities of the form: 

7 

= cTopo(ni,n2,n3,n4,n5,n6;n7) -F^ni(iiTopo(ni,- ■ ■,ni + l,- ■ -,^7) 

i=l 

1,7 

+ ^nieijTopo{ni, ■ ■ ■ , n^ + l, ■ ■ ■ , -1, ■ ■ ■ ; 77,7). (48) 

For a given set of the indices {nj}j=i 7 there are eight of such identities; in general, 
however, they are not all linearly independent. Each identity contains, in general, 
three kinds of amplitudes: 

1. the amplitude itself: Topo (ni, n2, na, 714, ns, ng; ^7); 

2. amplitudes with one of the indices increased by one, rij — >■ + 1 = 1, 7): 
Topo {rii, ■ ■ ■ , nj + 1, ■ ■ ■ ; n-i). These terms clearly originate from the derivation 
of the factor —niP~^'~^ . We pulled a factor Ui out of the coefficient 
di to point out that these terms are absent for rii = 0. In other words, it is not 
possible to generate an absent denominator by differentiation. This implies 
that the ibp identities for a given amplitude with topological number t involve 
only amplitudes with smaller or equal topological number, t' < t; 

3. amplitudes with one index increased by one, rii ^ rii + 1 and another index 
decreased by one, rij — >• rij — 1: Topo(r;,i, ■ ■ ■ , + ■ ■ ■ , n^ — 1, ■ ■ ■ ; 77,7). These 
terms originate from the derivation of the factor and the cancellation of 
a power of Pj by the invariants generated at the numerator. We pulled out 
a factor rij also out of the coefficients Cjj to point out similar properties as 
those discussed in 2). 

Other identities are obtained by considering discrete symmetries of the auxiliary 
diagram. The ladder topologies in (a), (b) and (c) of Fig. |3]for example have an 
up-down symmetry, resulting in the exchange of the fermion lines expressed by the 
relation: 

Topo (ni, ^2, ns, 77,4, n5, ng; n-j) = Topo {rii, ris, ^2, n4, uq, n^; n-j) . (49) 
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The massless crossed ladder also has an up-down symmetry, which interchanges also 
the boson lines. By solving the above identities, it is possible to express a general 
scalar amplitude as a linear combination of a finite set of Np basic integrals, called 
master integrals (Mi's): 



Nf 



Topo(ni, n2, n^, rig, rig; Uy) = ^ Ci(ni, na, ng, n^, n^, uq; uj) Fi. 



(50) 



i=l 



In more formal terms, we may say that the ibp identities induce such a strong linear 
dependence between all the infinite scalar amplitudes to project them into a finite 
dimensional linear space. The Fj's have the same role as the basis vectors in a linear 
space. The choice of the Mi's is arbitrary, as it is the choice of the basis vectors; 
only their number Np is fixed and equals the dimension of the space spanned by any 
choice of the -F/'s. Changing the set of Mi's in the decomposition is analogous 
to a change of basis. To give an explicit example, let us report the MI decomposition 
of the basic six-denominator amplitude of the ladder with a massive external boson 
in Fig. n(b), the example in the previous section: 



,2 7 \ f I 1 1 
Topo(l,l,l,l,l,l;0) = 2 --- + 6 [^ + - + 



X 



xTopo(0, 1,0, 1,0, 1;0) + 

1 



1 4 
^-- + 3 



1 



llli 

X 



1 — a; 



e J [l-xy 
[Topo(l,0,0,l,l,0;0) 
4 



-Topo(l, 0, 0, 0, 1, 1; 0)] -2 ( -- -+4) Topo(l, 1, 0, 1, 0, 1; 0) 



--2 I -Topo(l,0,0, 1,1,1; 0)- 
e I X 



2 9 



-9 



+ 



11 

e 



12 



1 



1 



X 



{i-xY 

Topo(l,0,0,l,l,l;-1) 



-2 ) 1topo(1, 1, 1, 0, 1, 1; 0)+ -Topo(l, 1, 0, 1, 1, 1; 0), (51) 
e x X 



where we have taken for simplicity a = 1. In this case Np = 8. We have chosen 
Mi's with a subset of the denominators with unitary indices: rii = 1, the remaining 
denominators having vanishing indices: rij = 0. For the topology represented by 
the amplitude Topo(l, 0, 0, 1, 1, 1; 0) this was ot sufficient as two master integrals 
are involved. We added the amplitude with the auxiliary denominator brought to 
the numerator, i.e. with index = —1: Topo(l, 0, 0, 1, 1, 1; — 1); we could have 
taken n2 = —1, = 0, n-i = or n2 = 0, = —1, = as well. Still another 
possibility would have been to consider amplitudes with one of the denominators 
squared, such as for example Topo(2, 0, 0, 1, 1, 1; 0). In general, the choice of the Mi's 
is dictated by practical considerations. In some cases, for example, it may be useful 
to require the absence of ultraviolet or infrared singularities from the Mi's ^Hl; in 
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other cases, the set of Mi's can be chosen in such a way that the corresponding 
system of differential equations become easier to solve (see next section). 

Let us now discuss the methods to solve the ibp identities to arrive to the results 
(j^Uj) . The ibp identities constitute a system of linear equations in which the un- 
knowns are the scalar integrals themselves. The oldest approach, developed in the 
original article on the ibp's, involves a symbolic solution of the identities, treated 
as recurrence equations on the indices [3 CHI- introduces operators raising or 
lowering one of the indices 

Topo(ni, ■ ■ • , • ■ ■ ; ny) = Topo(rai, ■ ■ ■ , ± 1, ■ ■ ■ ; riy). (52) 

Each ibp equation involves the identity operator 1 and the raising operators and 
i+j^. If one can combine the equations in such a way that an amplitude is written 
in terms of amplitudes all containing lowering operators, then reduction to simpler 
topologies has been achieved. The advantages of the symbolic approach are the 
elegance and that its implementation does not require in general a lot of CPU time. 
The main disadvantage is that a careful, case-by-case, inspection of the equations 
is required. A more recent approach is that of replacing numerical values for the 
indices and solve the resulting linear system OE]. Let us describe this method 
in practical terms. An initial linear system can be obtained by assigning to the 
indices the values = 0, 1 for i = 1, 6 and ny = 0. This is a set of 8 x 2^ = 512 
equations^. If these equations are not sufficient for a complete MI reduction, one 
may solve a larger system in which one of the indices can also take the value rij = — 1 
for z = 1, ■ ■ ■ , 7; this is a set of (3 x 2^ x 6 + 2'^) x 8 = 5632 equations. Let us stress 
that all the equations with one index = — 1 have to be considered simultaneously, 
because the amplitudes involved are coupled by the operators; for example: 

7+6" : Topo(l, 1, 1, 1, 1, 0; -1) ^ Topo(l, 1, 1, 1, 1, -1; 0). (53) 

One may also consider a system in which one of the indices takes the value rii = 2 
for i = 1,6; this amounts to3x2^x6x8 = 4608 equations. If these systems are not 
sufficiently large, one may consider equations with indices taking both the values 2 
and —1 and so on. The general idea is that, by going to larger systems, the number 
of equations grows faster than the number of amplitudes, till a "critical mass" of 
equations is reached, allowing the full reduction to MI's^" ^7]. The linear system 
can be solved with the method of elimination of the variables. In general, one has 
to devise a rule about which amplitudes in the system to solve first. Important 
quantities for the scalar integrals are the sum of the positive indices, 

7 

d = J2n,e{n,) (54) 

i=l 

^As already said, not all the equations in the system are linearly independent on each other and 
some of them are actually trivial (i.e. = 0), such as for example those with five indices equal to 
zero. 

^°One of us (U.A.) wishes to thank S. Laporta for a discussion on this point. 
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and minus the sum of the negative indices, 



7 



S = 



-rii d{-ni). 



(55) 



i=l 



A convenient criterion is to solve for the amphtudes with the greatest t, d and s first. 
The advantage of the numerical indices method is that it does not require a hand 
inspection of the equations — requiring patience and often a good mathematical 
intuition — and is therefore very general. The main disadvantage is that it may 
require a large CPU time and very long intermediate expressions may be generated 
because of the division by the coefficients of the unknowns amplitudes. 

Let us now discuss the ibp identities in the framework of the second reduction 
scheme discussed in the previous section, the one involving independent scalar prod- 
ucts at the numerator and loop momentum shifts. As with the auxiliary diagrams, 
the scalar integrals of Eq. (jmi) can be related by the ibp identities, which now read: 



where, as before, i = 1,2 and = ki,k2,Pi,P2- For a given set of indices {si,nj} 
of the input amplitude, we find eight identities involving integrals with up to an 
additional power on a denominator and an additional power of a scalar product at 
the numerator. Sub-topologies with t — 1, t — 2, ... coming from simplifications 
between scalar products and corresponding denominators also appear. As in the 
auxiliary diagram case, discrete symmetries of various topologies generate additional 
identities to be combined with the ibp identities 

Concerning our problem, we can reduce the calculation of all the scalar integrals 
related to the topologies shown in Figs. El and [7| — and therefore the calculation 
of the Feynman diagrams of Figs. Hand|21 — to that of the 22 Mi's in Fig. |H| Four 
of them are the product of two one-loop Mi's. The picture of the diagram alone 
denotes the scalar integral — only denominators — while the picture of the diagram 
accompanied by a scalar product denotes the scalar integral with that scalar product 
at the numerator. 

The algorithms explained in this section were implemented in a chain of programs 
written in the algebraic computer language FORM JHI and, in the case of the auxiliary- 
denominator scheme, also in Mathematica [201 • For the solution of the linear system 
of ibp identities, we used SOLVE This program uses the language C as a shell 
to manage a user-provided basic code written in FORM. External calls to Maple [221 
are done in order to simplify the coefficients coming in intermediate steps. 

3 Evaluation of the Mi's: the differential equa- 
tions method 

Let us consider a general amphtude Topo (ni, n2, ns, n^, n^, uq, n-j) in the framework 
of the auxihary diagram. As shown in the previous section, it can be decomposed. 




(56) 
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by means of the ibp identities, into a superposition of a set of primitive amplitudes 
called master integrals (Mi's). The latter cannot be computed with the ibp iden- 
tities, which, in general, only allow for a reduction of the number of independent 
amplitudes to be individually computed. The Mi's have been computed with a va- 
riety of techniques over the years: Feynman parameters, dispersion relations, small 
momentum expansions, large momentum expansions, just to mention the more com- 
mon ones. A technique developed over the last few years which is very efficient in 
our case is that of the differential equations in the external kinematical invariants 
[HI 13 Uni HH nil Uni mi ■ Let us then review the basics of this method. Once again, 
let us consider a specific example: the computation of the Mi's in Fig. |H1 (k) in the 
auxiliary diagram scheme. Let us define: 

Fi(s,a,e) = Topo(l, 0,0, 1,1,1; 0), (57) 
F2(s,a,e) = Topo(l,0,0,l,l,l;-1), (58) 

The above two Mi's appear in the decomposition of the six-denominator amplitude 
presented in the previous section (Fig. [11(b)). They are functions of s and a, as well 
as of the space-time dimension D = 4 — 2e. Let us take the derivative of Fi with 
respect to s at fixed a and with respect to a at fixed s: 

d d d 

s— Fi(s, a, e) = p'^^Fi{s, a, e) = p^-^Fi{s, a, e) , (59) 

d 

a—Fi{s,a,e). (60) 
oa 

The derivatives with respect to the external momentum component Pi or ^ind 
with respect to the mass squared a can be taken inside the loop integral, because 
of the properties of Dimensional Regularization. The resulting amplitudes can be 
reduced to combinations of Mi's according to the methods described in the previous 
section, so that: 

d 

s—Fi{s,a,e) = "^djis, a, e) Fj{s, a, e), (61) 
j 

d 

a—Fi{s,a,e) = '^Kij{s,a,e) Fj{s,a,e). (62) 
i 

Note the double role played by the ibp identities in the calculation: they allow for 
the reduction of independent scalar amplitudes into combinations of master integrals 
and they also allow for the generation of the differential equations to be solved for 
the Mi's themselves. 

Eqs. ()61|) and ()62|) can be simplified changing variables from s and a to x = —s/a 
and a. This triangular change of variables implies that: 

d d 
x—Fi{x,a,e) = s—Fi{s,a,e), (63) 

d d- 

a—Fi{x,a,e) = s—Fi{s,a,e) + a—Fi{s,a,e) = ^Fi{x,a,e). (64) 
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The second equation holds because of the dimensional relation 



(65) 



where di is the mass dimension of the MI: di/2 = D — ^^l,^^ n^^ . In the new variables, 
the scale a represents an over-all scale of the amplitudes and its evolution equation 
is consequently trivial. We therefore concentrate on the x-evolution equation only 
and we set a = 1 from now on in this section: 

Nf 



d_ 
dx 



Fi{x,e) = ^Aij{x,e)Fj{x,e). 



(66) 



In our case, by performing the above steps, one obtains the following linear system 
of coupled ordinary differential equations with variable coefficients: 

dFi 



dx 



dx 



{x,e) 



-Fi{x,e) 

X 



-a 



(2 - 3e) 
1 



1 



X 



[x,e 



X 



-a-e) 



1 — X 

2 - 3e 
1 



1 — X 
[F,{e) + F,ix,e)] 
1 



(67) 



1 

-+ . 
X 1 — X 



1 — X 

Fsie) - (1 - e) 



1 — a; 



F4{x,e) 



(68) 



where we defined: 



Fsie) = Topo(l, 0,0, 1,1,0; 0), (69) 
Fi{x,e) = Topo(l,0,0,0,l,l;0). (70) 

F3 is represented in Fig. |H1 (u); it is a sunrise diagram with one massive line, 
evaluated at a light-cone external momentum, equivalent to the null momentum. F3 
is then effectively a vacuum amplitude, so we dropped the dependence on x. It is 
obtained from Fig. |S1 (k) by shrinking one of the internal lines having an endpoint 
in the hard vertex. F4 is represented in Fig. |H1 (v) and is the product of a massless 
bubble and a tadpole; it is obtained from Fig. |H1 (k) by shrinking the massless line in 
the bubble. Note that there are no other non vanishing sub-topologies. We consider 
the sub-topologies F3 and F4 as known amplitudes'^, so that the system of Eqs. (jU7|) 
and (|UH|) is schematically written as: 



np 



^jyx,e) Fj{x,e) + Vli{x; e) 



(71) 



where the Qi{x, e) are known functions constituting the non- homogeneous part of the 
system. Let us note that the coefficients of the homogeneous terms have a regular 
expansion around e = 0, involving, in this case, only two terms: 



Aij(x, e) 



A^\x)+eAl/[x} 



(72) 

^^Their expressions are given respectively in Eqs. (|121H12^ (where we put F3 — F'^''^) ans Eqs. 
()11^8Hl^lj|l (where = F^^)). 
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A great simplification is offered by the fact that the system is triangular in our basis, 
because the equation for dF2/dx depends only on F2 and not on Fi. ^21(2;, e) = 0. 
Let us then consider this equation first. We substitute the explicit expression for 
the expansion of VI2 in powers of e and the symbolic expansion for F2, up to the 
desired order n: 

n 

n2{x,e) = J2''^2\^) + 0{e-+') , (73) 
i=-2 

n 

F2{x, e) = J] e^F^'\x) + O (e'^+^) . (74) 

The double pole — the leading singularity — obeys a differential equation of the 
form 

^(x) = Af^{x)Ft'\x) + n(-^\x) (75) 

where 

4?(x) = -i--^. (76) 
X 1 — X 

The associated homogeneous equation, 

duj' 



x)^Af^{x)u;2ix) (77) 



dx 

is solved by separation of variables to give: 

\2 



u!2{x) = exp 



A^^{t)dt 



X) 



X 



(78) 



with an overall constant omitted in the last member. The double pole coefficient is 
then computed with the method of of the variation of the constants: 

Ft^\x) = u;2{x) J' K2{t)n^-^\t)dt , (79) 

where the primitive involves an arbitrary constant to be fixed imposing initial con- 
ditions, and where the kernel K2 is given by 

Higher orders are computed through the recursive relation in e: 

= A^^ix)F^\x) + AS^ix)Ft'\x) + 0«(x). (81) 

The solution of the above equation is again obtained with the method of variation 
of constants: 

F^\x) = U2ix) J' K2{t)Q^'^{t)dt , (82) 
22 



where 

fi«(t) = n»(t) + Ag(t)F^^)(t) (83) 

is the inhomegeneous term redefined to include the amphtude of previous order in 
the e expansion. Note that the kernel K2 is the same for each order in e. 

The full evaluation of F2 requires, as we already said, an initial condition on the 
X evolution. The latter can often be derived by studying the behaviour of the MI 
close to a threshold or a pseudothreshold. In our case, let us consider the point 
X = 1: it is a pseudothreshold as it corresponds to s = —w? < and therefore is 
not a singular point for the amplitude, so that: 

.dF2(x,e) , 
lim 1 - x) — \^ = 0. 84) 
ax 

By imposing the above condition to the differential equation, an initial value for the 
amplitude in x = 1 is obtained: 

(2 - 3e)F2(l; e) + (1 - e)F3(e) + (1 - 6)^4(1; e) = 0. (85) 

Once F2 has been explicitly computed, one replaces the result in the first equation 
and solves it for F\ in the same way. 

In general, however, the system is not triangular and the above method is not 
directly applicable. One can try a transformation to a new set of Mi's, 

Fi, F2 ^ F[, F^, (86) 

in order to have the system in a triangular form. A triangularization to order e is 
actually sufficient, i.e. one for which J^^l = but A^^^ 7^ 

The Mi's can be expressed in terms of the one-dimensional harmonic polyloga- 
rithms (HPL's), the latter having the following recursive integral definition 

H{l;x) = -log(l-a;) (87) 
H(0,x) = log(x) (88) 
H{-l;x) = log(l + x) (89) 

and 

H{a,w]x)= f{a;y)H{w;y)dy, (90) 
Jo 

where the basic weight functions are defined as: 

/(i;^) = ^^^^ 

f{0;x) = - (92) 



^^In more complicated cases, it may happen that a triangular form of the system cannot be 
reached, and the coupled system oi np equations must be solved. This is equivalent to the solution 
of a differential equation of order np for anyone of the Mi's, for which no general algorithm exists. 
In our problem, the number of Mi's for a given topology is at most np = 2, corresponding to 
a second-order differential equation. Moreover, it turned out that the solution of the associated 
homogeneous system was rather simple, so that the complete solution could be found with the 
method of the variation of arbitrary constants. 
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For a complete discussion of HPL's and their properties we refer to [2]; for their 
numerical evaluation see 0. 

As an "empirical" rule, we found that it was sufficient to stop the e expansion 
beyond the order involving weight w = 4 HPL's. 

The massless Mi's cannot be computed with the differential equation method and 
a more traditional method as that of Feynman parameters or dispersion relations 
has to be used. As discussed in the introduction, they are indeed of the form: 

(massless amp.) = a^"""^ ( ^ j G(e)x'-"'^-2^ (94) 

Since the x-dependence is factorized, the x evolution gives only a trivial, dimensional, 
information, analogous to the a evolution. The non-trivial quantity to compute is 
G{e). In other words, all the information is contained in the initial data — for 
example the amplitude at x = 1 — that can be extracted only by direct integration. 

For a graphical summary of the whole method discussed in sections |21 and |2l see 
the flowchart in Fig. |H1 



4 Results for the master integrals 

In this section we present the results of our computation of the Mi's involving up to 
five denominators included, which constitute a necessary input for the calculation 
of the two-loop vertex diagrams in Figs. ^ and El They are expanded in a Laurent 
series in 

e = 2 - D/2, (95) 

up to the required order in e. The coefficients of the series are expressed in terms 
of HPL's [21IS1 of the variable x, defined as : 

X = — , (96) 
a 

where q = pi + p2, and^^ a = m? . We denote by /x the mass scale of the Dimen- 
sional Regularization (DR) — the so-called unit of mass. The explicit values of the 
denominators are given in appendix El We work in Minkowski space and we 
normalize the loop measure as: 

This definition makes the expression of the one-loop tadpole — the simplest of all 
loop diagrams — particularly simple (see appendix lB|) . 

^•^Note that with our definition of the scalar product s = —q^. 
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Scalar Amplitudes 



Independent 
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Master Integrals 



System of Diff. Eqs. 



Projection on invariant form factors 



Auxiliary diagram or sinifts 



IBP's and general symmetry relations 



Generation of Diff. Eqs. (+ IBP's) 



Solution of the System in Laurent series of e 



Figure 9: Flowchart of the method used for the reduction to the Mi's and their 
evaluation. 



25 



4.1 Topology t = 3 
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The above amphtude is easily computed with Feynman parameters or with an iter- 
ation of the general formula for a massless bubble given in 0. 
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where: 
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~2 ' 

_3 
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-^-C(2), 
-Y-3C(2) + C(3), 
-^-7C(2)-^ + 3C(3). 



The above amplitude is easily computed with Feynman parameters. 
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.(5) 



(5) 



.(5) 



= -3 + H{0,x), (129) 

= -7 + C{2) + 3H(0,x)-H(0,0,x), (130) 

= -15 + 3C(2) + 2C(3) + (7 - a2))H{0, x) - 3H{0, 0, x) 

+H{0,0,0,x), (131) 

= -31 + 7C(2) + ^C'(2) + 6C(3) + (15 - 3C(2) - 2C(3))i/(0, x) 

-(7 - C{2))H{Q, 0, x) + 3i^(0, 0, 0, x) - H{0, 0, 0, 0, x) . (132) 



4.2 Topology t = 4 
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1 + 



[4i7(-l,a;) -2if(-l,-l,x) 



-i^(-l,0,a;) + i/(0,-l,a;)], 
65 

- + 3C(2)-C(3)-7i/(0,x)-12 



(135) 
(136) 

(137) 



1 + 



H{-l,x) + 2H{0,0,x) 



+ 



x 



[4:H{-1, 0, x) - 4:H{0, -1, x) + 8i/(-l, -1, x) 

-AH{-1, -1, -l,x) - 2H{-1, -l,0,a;) + 2i7(-l,0, -l,x) 
-2H{-1, 0, 0, x) + 2i/(0, -1, -1, x) + i/(0, -1, 0, x) 
-i/(0,0, -1,,t)], 

^ + 5C(2) + - 9C(3) - (33 - 2C(2))i/(0, x) 



(138) 



-2(16 -3C(3)) 
1 



5 

1 

1 + - 

X 



H{-1, x) + 14if(0, 0, x) - 4if(0, 0, 0, x) 



+ 



1 + 



X. 



[2(6 - C(2))//(-l, 0, x) - 12H{0, -1, x) + 24//(-l, -1, x) 

-8//(-l, -1, 0, x) + 8H{-1, 0, -1, x) - SH{-1, 0, 0, x) 
+8H{0, -1, -1, x) + 4//(0, -1, 0, x) - 4//(0, 0, -1, x) 
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-16i/(-L -L -1, ,t) + 8i/(-l, -1, -1, -1, ,t) 
+4if(-l, -1, -1, 0,x) - 4:H{-1, -L 0, -1, x) 
+AH{-l,-l,0,0,x) -AH{-1, 0,-1,-1, x) -2Hi-l, 0,-1,0, x) 
+2H{-1, 0, 0, -l,x) + AH{-1, 0, 0, 0, x) - AH{0, -1,-1, -1, x) 
-2H{0, -1, -L 0, ,x) + 2H{0, -1, 0, -l,x)- 2H{0, -1, 0, 0, x) 
+2H{0, 0, -1, -1, x) + H{0, 0, -1, 0, x) - H{0, 0, 0, -1, x)] . (139) 



^ ^ i=-2 



where: 



f!5 = 1, (142) 

F^j = A-2H(0,x), (143) 

FjS^^ = 12-2C(2)-8i/(0,x)+4i/(0,0,x), (144) 
Fi^') = 32-8C(2)-4C(3)-4[6-C(2)]i/(0,x) + 16if(0,0,a;) 

-8i/(0,0,0,x) , (145) 

= 80 - 24C(2) - - 16C(3) - 8[8 - 2C(2) - C(3)]if(0, x) 

5 

+8[6 - ((2)1^^(0, 0, x) - 32H{0, 0, 0, x) + 16H{0, 0, 0, 0, x) . (146) 



= (^)"e^'';'''+o(^'). 

^ ^ i=-2 



Pop 



9-'^12 



where: 



.(8) 



1, 

4-/^(0,2;) 



1 + 



x 



H{-l,x), 
1 



Fo^°^ = 12 - C(2) - 4//(0, x) - 4 1 + - H(-l, x) + H(0, 0, x) 



X 



+ 



1 

1 + - 

X 



[H{-l,0,x) + 2H{-l,l,x)], 



(147) 
(148) 

(149) 
(150) 

(151) 
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X 



X 



32 - 4C(2) - 2C(3) - [12 - C(2)]//(0, x) - [12 - C(2)] [l + 

[4H{-l,0,x) + 8H{-l,-l,x)] 
[H(-l, 0, 0, x) + 2H(-1, -1, 0, x) 



xH{-l,x)+AH{0,0,x 

1 



1 

1 + - 

X. 



-H(0,0,0,x) 
+4Hi-l,-l,-l,x)], 

80 - 12C(2) - 



x. 



(152) 



10 



-8C(3)-[32-4C(2)-2C(3)] H{0,x) 



+ (l + ^)//(-l,x)] + [12 - C(2)]{//(0,0,x) + (^1 + 1) X 
x[2i/(-l,0,a;) + 2i/(-l,-l,a;)]} - 4i/(0, 0, 0, 



a; 



[8ff(-l, -1, 0, x) +4//(-l, 0, 0, x) + 16//(-l, -1, -1, x) 

//(-I, 0, 0, 0, x) - 4:H{-1, -1, -1, 0, x) - 2H{-1, -1, 0, 0, x) 
SH{-1, -1, -1, -1, x)] + H{0, 0, 0, 0, x) . (153) 





(P2 • h 



^2(4-D) {dDk,}{d''k2} 



^ ^ i=-2 

(^)"^,.f;M+0(e^), 



12 



(154) 
(155) 
(156) 
(157) 



1 

2' 

l-H{0,x), 

19 r 1" 

--2C(2)-5//(0,x) + //(0,0,x)- 1-- i/(l,0,x) 
+ -i/(0,l,0,x), 

X 

fi5 

y - 10C(2) - C(3) - [19 - C(2)]//(0, x) - 3C(2) 



(158) 
(159) 



(160) 



1 

1 - - 

X 



H{l,x) 



+5H{0, 0, x) + ^^^H{0, 1, x) - 5 



X 



1 

1 - - 



H{l,0,x) 



31 



-H{0, 0, 0, x) + -H{0, 1, 0, x) - 



X 



1 - 



X 



[^(0,1, 0,x) 



-H{1, 0, 0, x) + 3H{1, 1, 0, x)] + -[H{0, 0, 1, 0, x) 



X 



-H{0, 1, 0, 0, x) + 3H{0, 1, 1, 0, x)] , 



X 

16' 

^[9-AH(0,x)], 

^ + ^[63 - 16C(2)] - ^[2 + 3x]H{0, x) + ^i/(0, 0, x) 



1 

x 

X 



i/(l,0,a;) + -iff(0,l,0,a;), 
4a; 



57 405a; 9, a; 1 



-4C(2)x]i/(0,a:)- 



3C(2) 



(161) 
(162) 

(163) 

(164) 



1 r 



11 



X - - H{1, a;) + — [2 + 3a;]ii"(0, 0, x) 
X -J j-Ky 

H{1, 0, x) + ^^^(0, 1, x) - ^H{Q, 0, 0, x) 



[i/(0, l,0,a;) -i/(l,0,0,a;) 



i« 2 + 9a;- 
16 L X 

+^//(0,1.0.x)-i^x-^ 

+3iJ(l, 1, 0. a;)] + —[H(0, 0, 1, 0, x) - H(0, 1, 0, 0, a;) 

4a; 

+3/^(0, 1,1,0, a;)]. 



(165) 




(166) 
(167) 



1 

2' 
19 

T 

65 



-5H{0,x) + 2H{0,0,x), 

- 4C(3) - 19if(0, x) + 10i/(0, 0, a;) - AH{0, 0, 0, a;) , 



91 1 19 

— - -C'(2) - 20C(3) - [65 - 8C(3)]//(0, x) + 38H{0, 0, x) 
-20//(0, 0, 0, x) + 8H(0, 0, 0, 0, a:) . 



(168) 
(169) 
(170) 

(171) 

(172) 
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As discussed in section 3, the above amplitude, being massless, cannot be calculated 
by means of the differential equations technique. It is easily computed with Feynman 
parameters. 



(173) 
(174) 




where: 



.(12) 



^(12) 



.(12) 



^(12) ^ 



1 

2' 
3 

2' 



1-1 

X 



H{l,Q,x), 



(175) 
(176) 
(177) 



■- + 5C(2)-C(3) + 7i/(0,a;) + 2C(2) 



1-1 

X 



H{l,x) 



-2H{0,0,x) + 



1 

1 - - 

X 



[AH{l,0,x) + H{0,l,0,x) 



.(12) 



-2i/(l,0,0,a;) + 2i7(l,l,0,a;)], 

+ 19C(2) + ^ + C(3) + [33 - 2C(2)]i/(0, x) + 4[C(3) 



(178) 



+2C(2)] 



+2[6 - C(2)] 



X 



H{l,x) - UH{0,0,x) + 2C{2) 



X 



H{0,l,x) 



1 



X 



+AH{0,0,0,x) + 



1 - 



H{l,0,x)+4C{2) 
1 



1 - 



1 



X 



H{l,l,x) 



X 



[AH{0,l,0,x)-8H{l,0,0,x) 



+8H{1, 1, 0, x) + H{0, 0, 1, 0, x) - 2H{0, 1, 0, 0, x) 
+2H{0, 1, 1, 0, x) + AH{1, 0, 0, 0, x) + 2H{1, 0, 1, 0, x) 
-4:H{1, 1, 0, 0, x) + AH{1, 1, 1, 0, x)] . 



(179) 



4.3 Topology t = 5 




= (^)"e^"^<""+ow. 

^ ' i=—l 



8-^^13 



(180) 

(181) 
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where: 



aF, 



(13) 



-[//(0,0,-l,x) + //(0,l,0,x)], 

X 

-[3C(2)ii'(0, 1, x) - 4H{0, 0, -1, -1, x) + H{0, 0, 0, -1, x) 
+H{0, 0, 1, 0, x) - H{0, 1, 0, 0, x) + 3H{0, 1, 1, 0, x)] . 




where: 



aF^'^^ = -[C(2)i/(0,l,x) + //(0,l,0,0,x) + //(0,l,l,0,x)] 

.37 




where: 



(15) 



aFr 



aF 



(16) 



aFr 



(16) 



(16) 



ifF^^ + Oie), 
a J 




{h ■ k2) 



= (^)"E^'^;"°' + f(.=), 

^ ^ j=-2 



A:2 



-[2C(2)i/(0 -1, x) -i/(0 -1, 0-1, x)+i/(0 -1, 0, 0, x)] , 
1 

2 ' 

2-i^(0,x), 

7-3C(2)-^//(0,x)-(2-C(2)) 



+^/7(0,0,a:) + i 



3-- 



1 



X 



1 + -]h{-1,x) 
H{0-l,x) + l\l + -][H{-l,0,0,x) 



-H{-l,0,-l,x)], 

24 - ^C(2) + ^ - ^[7 - 3C(2)]//(0, x) - (l3 - 2C(2) 
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+|C(3)) 



1 + 



X 



H{-l,x)- 



C(2) 



1 + 



H{-l,0,x) 



+5i/(0,-l,x)+ C(2)- 



1 + 



H{0,-l,x) 



+2(4 -C(2)) 



1 

1 + - 

x. 



H{-l,-l,x) + -HiO,0,0,x) 



-8H{0-1-1, x) + 2H{0, 0,-1, x) + - 



X 



[4H{0~l-l,x) 

-H{0,0,-l,x) - 2H{-l,0,-l,x) +2H{-l,0,0,x) 
+2H{-1-1, 0,-1, ,T)-2iJ(-l,-l, 0, 0, x)+AH{-l, 0,-1,-1, x) 
-H{-1, 0, 0, -l,x)- 3H{-1, 0, 0, 0, x) - H{Q, -1, 0, -1, x) 
+H{0,-l,0,0,x)]. (195) 

The above topology is the only 5-denominator one having 2 Mi's. 



= (^)"E^'';""+0W. (197) 

^ ' j=— 1 




where: 



(17) 



-H{Q,l,0,x), 

X 

-[2C{2)H{0, 1, x) + H{0, 0, 1, 0, x) - 2H{0, 1, 0, 0, x) 



+2H{0,l,l,0,x)]. 



(198) 



(199) 




where: 



= (^)"e^*^;'"*+o(^=), 

^ ^ i=— 1 



1 



^ --H(i,o,x), 

X 

1 



(200) 
(201) 

(202) 



(18) 







aFi''^ = --[2(C(2)+C(3))//(l,x)+C(2)//(l,l,x) + (4-3C(2))//(l,0,x) 



— [C{'2)H{1, x) + 2H{l, 0, x)-3H{l, 0, 0, a;) + i/(l, 1, 0, x)], (203) 



-6/7(1, 0, 0, x) + 2i/(l, 1, 0, x) + 7i/(l, 0, 0, 0, x) 
-3H{l,l,0,0,x) + H{l,l,l,0,x)]. 



(204) 
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The above amplitude has a simple ultraviolet pole coming from the massless 
bubble. 



^ ^ i=0 



(205) 
(206) 



where: 



,(19) 2 



aF^'"" = -[C{2)H{l,x)-H{0,0,-l,x) + H{l,0,0,x) 



X 

-2H{l,0,-l,x)], (207) 
aF^ = -[(2C(2)-C(3))i/(l,x) + C(2)(if(l,0,a;) + iJ(l,l,a;)) 

-2H{0, 0, -1, x) - 4/7(1, 0, -1, x) + 2H{1, 0, 0, a;) 
+4:H{0, 0, -1, -1, a;) - H{0, 0, 0, -1, x) + 8H{1, 0, -1, -1, x) 
-3H{1, 0, 0, -1, x) - 3/7(1, 0, 0, 0, x) - 2i7(l, 1, 0, -1, x) 
+i7(l,l,0,0,x)]. (208) 

The finite part of the above amplitude is also given in 




;.^(^-^) / ^ (209) 

^ ^ i= — l 



where: 



aFlf) = --iy(l,0,x), (211) 

X 

aF(J'°^ = --[C(2)/f(l,x) + if(0,l,0,x) + 2i7(l,0,x) + i7(l,l,0,x) 

X 



-2H{l,0,0,x)], (212) 
aFf°) = --[(2C(2) + 3C(3))//(l,x) + C(2)i/(0,l,x) + 4f/(l,0,x) 

X 

+({2)H{1, 1, x) + H{0, 0, 1, 0, x) + 2i7(0, 1, 0, x) 

-2H{0, 1, 0, 0, x) + H{0, 1, 1, 0, x) - 4if(l, 0, 0, x) 

+4i7(l, 0, 0, 0, x) + H{1, 0, 1, 0, x) + 2i7(l, 1, 0, x) 

-2H{1, 1, 0, 0, x) + H{1, 1, 1, 0, x)] . (213) 

The above amplitude has a simple ultraviolet pole coming from the nested bubble. 



36 



5 Six- denominator diagrams 



In this section we present the results of all the six-denominator diagrams with up 
to one massive line entering the EW form factor. The ladder and vertex insertion 
diagrams — the planar topologies — are all reducible, while both the crossed lad- 
ders — the non-planar topologies — are Mi's, i.e. cannot be reduced to simpler 
topologies by using the ibp identities. Even though these diagrams have a different 
"diagrammatic status" , we present them together because they are the most impor- 
tant ones and may be used for reference. Series representations for the pole and 
finite parts of the ladder and crossed ladder diagrams have been given in j2SI 121] • 
The authors compute a large number of terms in a small momentum expansion of 
the diagrams and then fit the results to an assumed basis of functions and trascen- 
dent constants. In the above paper, formulas for the resummation of the series in 
terms of Nielsen polylogarithms [H |3] are also reported. By using them, we could 
check that the pole parts of our results agree with theirs. Finally, the leading loga- 
rithms of all the six denominator amplitudes have been computed in [22] by means 
of a pole approximation for the vector boson propagators; we are in agreement with 
them also. 

5.1 Planar topologies 



where: 




a^An 



/i 



2{4-D) 



{d^ki}{d''k2} 



Di 



/ 2\ 2e 
^ / i=-4 



1 

1 

~2^2 



H{0,x) 



x^ 



C(2)+if(0,0,x) 



5C(3) - 2C(2)i/(0, x) - 2H{0, 0, 0, x) 



-10C(3)if(0,a;)+4C(2)iJ(0,0,x) 



+4/7(0, 0, 0, 0, x) 



(214) 
(215) 

(216) 
(217) 
(218) 
(219) 

(220) 



The above diagram has also been computed in [26] using Feynman parameters 
and in [21] using dispersion relations; we computed it by means of the reduction 
to simpler (massless) amplitudes, the latter computed with Feynman parameters. 
We are in complete agreement with these references. The leading singularity of this 
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diagram is a pole in e of fourth order or, equivalently, a forth power in log(a;) in 
the finite part for x ^ 0. We have indeed a soft and a colhnear singularity for each 
loop, each one contributing by a simple pole to the amplitude. 




1 



2\ 26 



Bo + O (e) 



(221) 
(222) 



where: 



X' 



■[6C(3)i/(l, x) - 2C(2)i7(l, 1, x) + i7(l, 0, 1, 0, x) 
-2H{1, 1, 0, 0, x) - 2H{1, 1, 1, 0, x)] . 



(223) 



As expected on the basis of one's physical intuition, the above diagram is infrared 
finite. The massive vector boson cannot indeed propagate for large distances and 
the same is true for the photon, which is "trapped" inside the external loop. 




(-) E^'C. + fW. (225) 

^ ^ ^ i=-2 



where: 



a'C_2 = -^H{l,0,x), (226) 

a^C.i = -^[({2)H{l,x) + 4H{l,0,-l,x) + 3H{l,l,0,x) 

-3iJ(l,0,0,x)], (227) 

a'Co = -^[2C(3)i/(l,x) + 7C(2)i7(l,l,x)-3C(2)i/(l,0,x) 

-16i/(l, 0, -1, -1, x) + 6H{1, 0, 0, -1, x) + 7H{1, 0, 0, 0, x) 
+3H{1, 0, 1, 0, x) + AH{1, 1, 0, -1, x) - 5i/(l, 1, 0, 0, x) 
+9i/(l,l,l,0,x)]. (228) 

Unlike Eq. (j221|l , in the above diagram the photon is not trapped by the massive 
line and can propagate to asymptotic times. As a consequence, we have a double 
infrared pole, having a similar dynamical origin to the one in the massless one-loop 
vertex. 
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where: 



X {1 + x) 
1 1 



2(4-D) f {d^k^jid^'h} 

/ 2 \ 2e 

(^) i:^^E.^Oie). 



[C(2) + if(0,-l,x)], 
[C(3)-2C(2)(if(0,x) + /f(-l,x)) 



(229) 
(230) 

(231) 



X (1 + s) . 

-2H{-1, 0, -1, x) - 2i7(0, -1, 0, x) + H{0, 0, -1, x)] , (232) 

1 1 ir e{2) 



X (1 + x). 



10 



-C(3)[i/(0,x)+iJ(-l,x)] 



^(2) [2i/(0, 0, x) + 2i^(-l,-l, x) + 2/7(-l, 0, x) -3if(0,-l, x)] 
-4/7(0, -1, -1, -1, x) + 2/7(0, -1, 0, -1, x) + /7(0, -1, 0, 0, x) 

+ ^H{0, 0, 0, -1, x) - H{0, 0, -1, 0, x) + 2if(0, -1, -1, 0, x) 
+2i7(-l, 0, -1, 0, x) - //(-I, 0, 0, -1, x) 
+2/7(-l,-l,0,-l,x)}. 



(233) 




2{4-D) [{d^kAid^'k^} 



where: 



2x2 



[l + if(0,x)]. 



x^ 



:[l-C(2) + if(0,x) + i/(0,0,x)], 



a'Go = ^[2-2C(2) + 3C(3) + 2(l -C(2))i/(0,x) + 2i7(0,0,x) 



x^ 



+2i7(0,0,0,x)] 



(234) 
(235) 

(236) 
(237) 
(238) 

(239) 



The above diagram has also been computed in [26] using Feynman parameters 
and in j2ZI using dispersion relations; we have computed it by means of the reduction 
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to simpler (massless) amplitudes, the latter computed with Feynman parameters in 
section 3. We are in complete agreement with the previous results. The massless 
vertex correction has at most a triple pole in e and is therefore less singular than the 
massless ladder or crossed ladder (see later). As well known, the vertex correction 
is indeed neglected in the double logarithmic analysis of QED form factors. 




/ 2 \ 2e 
^ ^ i=-3 



where: 



a^I-3 



3 

1 

2^ 



[{l + x)H{-l,x)-xH{0,x)], 



:[C(2)x - 2xH{0, -l,x) + 2((1 + x)H{-l, -1, x)] , 



(240) 
(241) 

(242) 
(243) 
(244) 



X- 

^ |C(2) (1 + x)Hi-l, x) - mxH{0, x) - xHiO, 0, 0, x) 
+2a;[4i7(0,-l,-l,x) -i/(0,0,-l,x)] + (1 + a;)[if(-l, 0, 0, x) 
+2i7(-l,0,-l,a;) -8ii'(-l,-l,-l,a;)]} . (245) 




where: 



= .^'--/{<^"M{<^''M^J^^^^ (246) 



/ 2 \ 2e 
^ ^ j=-2 



(247) 



= (248) 

= ^[C(3)-2C(2)(i7(0,a;)+i/(l,a;))-2i/(0,l,0,a;) 

-2//(l, l,0,x)] , (249) 

aVo = -i{^ + C(3)[i/(0,x) + i/(l,x)]-C(2)[2i/(0,0,x) 

-3i/(0, 1, x) + 2i/(l, 0, x) - 3i/(l, 1, x) - //(O, 0, 1, 0, x) 
+2//(0, 1, 0, 0, x) - 3//(0, 1, 1, 0, x) - H{1, 0, 1, 0, x) 

+2i/(l,l,0,0,x) -3//(l, 1,1,0, x)}. (250) 
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/ 2\2c 



a'K. 




3 
3 

^[3 - 2C(2) - 2H{0, x) + 4i/(0, 0, x) 



[l-2HiO,x)], 



/ 2 \ 2e 



13 



L_2 = 



2a; 
1 

4a; 



a; 
17 

a; 



+ (i-i)ff(o..)-(i + ^ + ^)ff(-i..) 

+ 2C(2)(i-1) + (3-5)h(0,.)-3(i + ; 
+4)ff(-l.a;) - 2(1 + i)i?(0,0,i) -6ff(0,-l,x) 

+8(l + ^ + ^)^(-l'-l'^)- 




P J I 1/1 'ix>2VlV^V'! 



12 



1 

2^ ' 



41 



a^Mo 



--(7-2C(2)) + 
+2//(l,0,x)]. 



rl 
- + 



X {1 — x) 



12C(2) + ff(0,i)-2H (0,0,1) 




/ 2\ 2e 



where: 



-4(l + i^(0,0,x)) 




where: 



1 

1 



[l-^(0,x)], 



(264) 

(265) 

(266) 
(267) 

(268) 
(269) 
(270) 

(271) 
(272) 

(273) 
(274) 
(275) 



2x 

-[1 - H(0, x) + H(0, 0, x) - H(l, 0, x)] , 

X 

-[2 - 2C(3) - 2H{0, x) + C(2)i/(1, x) + 2H{0, 0, x) 

X 

+2H{1, 0, x) - 2H{0, 0, 0, x) + H{0, 1, 0, x) - 2H{1, 0, 0, x) 
+H{l,l,0,x)]. (276) 



(277) 
(278) 
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where: 



aP_2 = ^, (279) 

aP„i = -^H(0,x), (280) 

aPo = ^[l-C(2) + if(0,0,x)]. (281) 
x"^ 

The above amphtude is actually a 5 denominator one; we include it here just for 
practical convenience. 

5.2 Crossed ladders 

The massless crossed ladder is a MI, i.e. it cannot be reduced to simpler (massless) 
topologies by means of ibp identites. Its direct computation is rather involved so we 
report the result as taken from ref . |2S] • We checked Gonsalves result for the crossed 
ladder with Eq. (3.4) in [2Z]. If we include in the latter equation the dimensional 
factor a;"^"^*^ (in our notation), a complete agreement between the two computations 
is found^^. 



where: 




2(4-D) 



Dr 



1 



\ / j=_4 



X^ 



X^ 



;H{0,X), 

:[7C(2)-4iJ(0,0,x)] 



-^[27C(3) - 14C(2)if(0, x) + 8H{0, 0, 0, x)] 



x^ 



x^ 



57C^(2) 



54C(3)i/(0, x) + 28C(2)i/(0, 0, x) 



-16/7(0,0,0,0,3;) 



(282) 
(283) 



(284) 
(285) 
(286) 
(287) 

(288) 



Note that the above diagram is the only one having a fourth order (infrared) pole 
in e together, as we have seen, with the planar ladder. 

^"^Note that the definition of e in the above cited Van Neerven paper is different from ours: 
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1 



/ 2\2e 



where: 



a^Ro 



1 



1 



X (1 — x) 

1 1 
- + 



[C(2)-2i/(0,-l,a:)], 

X ' (i_^)j[C(3)-2C(2)(//(0,x)-2//(l,x)) 

+8H{Q, -1, -1, - 4/7(0, 0, -1, a;) + 2/7(0, 1, 0, a;) 
-8i/(l,0,-l,x)], 
1 I 1 If C^(2) 



(289) 
(290) 

(291) 



(292) 



-C(3)[i/(0,x)-2i/(l,a:) 



X (l-x)Jl 10 
+C(2)[2//(0, 0, x) + 8H(1, 1, x) - 4/7(1, 0, x) - 6//(0, -1, x) 
+7H{0, 1, x)] - 16//(0, -1, -1, -1, x) + SH{0, -1, 0, -1, x) 
-2iJ(0, -1, 0, 0, x) + 8H{0, 0, -1, -l,x)- 2H{0, 0, 0, -1, x) 
+/J(0, 0, 1, 0, x) - 8/J(0, 1, 0, -1, x) - H{0, 1, 0, 0, x) 
+3H{0, 1, 1, 0, x) + 16/7(1, 0, -1, -1, x) - 8/7(1, 0, 0, -1, x) 

+4//(l,0,l,0,x) - 16//(l,l,0,-l,x)}. (293) 

The differential equation for the above amphtude is a rather lengthy expression, as 
it involves fourteen subtopologies. 



6 Large momentum expansion, \s\ ^ m? 

In this section we give the asymptotic expansions of the six denominator scalar 
integrals, i.e. the expansion in powers of 1/x and L = log(a;) up to the order l/x"^ 
included. As already noted in the introduction, these expansions are relevant for 
the study of the structure of the infrared logarithms coming from multiple emission. 
The exact expression of the massless diagrams involves only an overall power of x 
and powers of L, so their expansion coincides with the former and is not duplicated. 
The coefficients not explicitly given are zero. 



where: 
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-10 - 2C(2) + 4C(3) - 4[1 + C(2)]L - - , 

2 D 



L2 



+ 2C(2) 



8 12 



(296) 
(297) 




2\ 2e 4 
/ 1=2 i=-2 



where: 



(29J 



a'Cz! 



2C(2) + ^L^ 



1 



3C(3)-C(2)L-|L^ 

^^-llC(3)L + 2C(2)L^ + iL^ 
10 

-1-L, 
C(2)-L + L2, 

-26 + 3C(2) + 11C(3) - 2[3 + 2C(2)]L + 3^^ - 
4 2 ' 

4 2 4 2' 
_117 ^ 17C(2) , 11C(3) 



4 



(299) 
(300) 

(301) 

(302) 
(303) 

(304) 
(305) 
(306) 



13 

T 



^ + 2C(2) L + ^L"- , (307) 



The expansion above is in agreement with that obtained in [28j by using a kinematic 
region decomposition of the integrand. 
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n 



i=2 



i=-2 



(308) 



where: 



a^EZo 



a'Ezl 



a 



a'E'. 



c(2) + ^L^ 

5C(3) 3C(2)^ 1^3^ 
4^-8a3)L+kML^ + ^L^ 



48" 



_1_^(2)-1l2, 

17 7C^(2) 
-- + 2C(2) - - 5C(3) - 4L + 3C(2)L + 8C(3)L 



(309) 
(310) 
(311) 
(312) 
(313) 
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^2 



1 



48 



5C(3) 3C(2) 
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where: 



where: 



2 4 
7C(2)^2_3 
4 



+ 



2 4 

29_7 C(2) , 7C'(2) , 15C(3) , 11 

32 

1 

2 



2 

7 

48' 



9C(2) 



L - 8C(3)L 




\ / 7=1 i=-3 V / 



_3 

~4' 

= -C(2), 
= 4C(3), 

^ ~2~2^' 
= 2 + 2L + L% 



a-'lzl 



a'Jzl 
a'JZl 



-5 + 3C(2) - 4C(3) - 5L + 3C(2)L - ^L^ ^ 
1 

"4' 



1 

2 

21 



+ L, 
I 3C(2) 



8 
1 

12' 

9 3 ' 
211 



+ -L--L\ 
4 4 



C(2) 47 



216 




2 36 



— L+— L^ 
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a ) ^ 

/ j=l i=-2 



+ o{- 



C(2) 



C(3) 



(314) 
(315) 
(316) 

(317) 
(318) 

(319) 

(320) 
(321) 

(322) 

(323) 

(324) 

(325) 

(326) 

(327) 

(328) 

(329) 

(330) 

(331) 

(332) 
(333) 
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a^JZl 



^2 T-1 



2 tO 



a J 



6C'(2) 



1 



l + C(2) + L + -L^ 

-4 + C(2) - 4C(3) - 4L + C(2)L - - '^L' , 

5 02) 
8 2 4 
13 , 7C(2) 



9 



C(2) 



2C(3) + , 2 



L + -L' 



59 ^ C(2) 



108 ' 3 18 ' 6 
67 29C(2) 4C(3) 155 



648 



+ 




18 3 



L + -L\ 
6 ' 

+ 



^ C 2 r 37^2 2 „ 
108 3 36 9 



/ 7=1 i=-2 



+ 01 - 



(334) 
(335) 
(336) 
(337) 
(338) 
(339) 
(340) 

(341) 



where: 



where: 






i 

~2' 


(342) 


^2 rO 

a^Lzl = 
a^Lzl = 
a'L% = 


-C(2j , 

2, 

2-3L, 

5 - 2C(2) -3L + 3L2 , 


(343) 
(344) 
(345) 
(346) 


a'Lzl = 


4 2 ' 


(347) 


a^L°3 = 


-2 + ^L + L^ 
4 


(348) 


a'Lzl = 


1 

~6' 


(349) 




1 2^ 
6 + 3^' 


(350) 




2e 4 / 1 \ 


(351) 




1 

~2' 


(352) 


a^Ml^ = 


1 

~2' 


(353) 




4-C(2), 
47 


(354) 



a^Mzl 
a^Mzl 



= -L, 

= 2C{2)-L + 2L\ 

= -L, 

= 2 + 2C(2) - 3L + , 

= -L, 

= ^ + 2C(2)-4L + 2L^ 




2\ 2c 4 



1 ° 



i=-3 



a^Ozl 

a^Ozl 
a'Ozl 

a'Ozl 



a^Ozl 




1 

4' 

2 2 ' 

1- 2C(2)-L, 

2- 4C(2)-C(3)-2L-C(2)L, 

-l-c(2)-L-^L^ 



1 1 r 

4 + 2^' 



C(2) 



+ -L - -L" , 



1 + 

36 3 6 2 



a) ^ ^ ^ ^ ( x5 

/ j=2 i=-1 \ 



_ki^3C(2)L-5L3. 



-1 + 



10 

C(2) 



2 2 



8 + 2C(2) - + 9L + 3C(2)L + 2L^ - ^L^ , 
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34 + 15C(2) - ^C'(2) - 42C(3) - 2L + 12C(2)L + 36C(3)L 



■13L 



10 

y 



13 

24' 



"4 + ^ + 2 ' 



Y + 3C(2) - ^ + f L + 3C(2)L + 3L^ - ^L^ , 
f + T-^-63C(3).^L.18C(2)..36C(3)L 



(378) 
(379) 
(380) 



_ - ^2 _ i^^2 _ ^ -^4 



39 

T 



13 

24' 



(381) 



7 Small momentum expansions, |s| <C 

In this section we present the small momentum expansions of the six-denomina- 
tor diagrams, i.e. expansions in powers of x and L = log(a;) up to first order in x 
included. As in the previous section, we do not consider the massless scalar integrals 
since their expansion is trivial. The coefficients not explicitly given are zero. 
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(382) 



where: 



a'B'i, = 6C(3), 
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1 



-3-C(2) + 3C(3) + 2L-^L2, 
-g-C(2) + 2C(3) + ^L-iL^ 



(383) 
(384) 

(385) 



where: 
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- C(2) - 


a'C\ - 
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- 3C(2) - 




1 

4 


-11, 

2 ' 




5 
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C(2) 
2 



2 6 



9 
4^ 



(387) 
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(389) 

(390) 

(391) 
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16 4 

9 3 ' 

31 C(2) 11 



49 ^ , 3C(2) ^ , 17 
"8 



- C(3) -^L+ + ^L^ - -^L' 



8 



7 

12" 



36 



6 



where: 



a Eq 



^2 pO 



a'E^' 
a^El 



where: 



4709 23C(2) 2C(3) 160 59 ^ 7 ^ 

"648" "~6 3~"^^ + ^^^^^+36^ " 18^ 



^ .7=-l j=-2 




02) 

2 
C(3) 



-C(2)L, 



-C(3)L + C(2)L^ 



53 3( 2 C3 5^ ^.^.^ 
16 2 2 4 ' ' 

-C(2)L^ 

333 29C(2) Ci2) 3C(3) 39, , , 



32 



10 



-^^^+c(2)L^ 




^ ^ i=— 1 i=— 3 



3 

"4' 

2 2 ' 



(392) 
(393) 
(394) 
(395) 

(396) 



_C!(2) 

10 
1_C(2) 
2 2 ' 

^-C(2)-^-L + C(2)L, 
I - 5C(2) + lc^(2) - C(3) - 3L + 2C(2)L + C(3)L (402) 



(397) 
(398) 
(399) 
(400) 
(401) 



= C(2), 



(403) 

(404) 
(405) 



(406) 
(407) 

(408) 

(409) 
(410) 
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„2 7—2 



l + C(2)-L-C(2)L+iL2-iL^ 
1 

~4' 

-1, 

1 C(2) 3^ 1 2 

1 

12 ' 
1 

2' 

^ _ M + 1 r _ 1 r 2 
216 6 36 12 ' 

^ ^ 7=-l i=-2 



,2 tO 



C(2) 
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C(3) 



C(2)L, 

-C(3)L + C(2)L^ 



10 
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9 - 5C(2) - C(3) - 5L + 2C(2)L + L^ 



C(2) 
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81 _ 11C(2) _ C(3) _ 27 3 
16 4 2 8 ^ 4 



2e 1 



j=-2 i=-2 



a^Lzl = 2, 
a^Lzl = 2-2L, 
a^L% = 2-2C(2) -2L + L2, 
1 

~2' 

4 2' 

2r0 7 C(2) 3^ 1 2 

^ 8 2 4 4 
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where: 



where: 




a^Mzl = 



1 

~6' 
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~2 ' 
1 

24' 
11 

48 ' 

2\ 2e 1 







j=-l i=-2 



2' 



2 
3 



-2 + 2C(2) + 3L - , 



a^M" = -- + 2C(2) + 4L-L2, 
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4 ' 

2 2 ' 

2-2C(3)-2L + L'-iL', 
-6 + C(2) + 5L - L^ 

4 2 4 2 



(433) 
(434) 
(435) 
(436) 

(437) 

(438) 

(439) 

(440) 

(441) 
(442) 
(443) 

(444) 
(445) 

(446) 
(447) 
(448) 

(449) 

(450) 
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^ ^ j=-l i=-2 
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where: 



a'RZl = (455) 

a'Rzl = ^-C(2)L, (456) 

a'R\ = -^-C(3)L + C(2)L^ (457) 

a'Ro' = -1 + ^, (458) 

aX"' = -4 + 2C(2) + ^ + L-C(2)L, (459) 

a^RO = _l4+5^(2)_i_Li+2C(3)+7L-4C(2)L-C(3)L--L^ (460) 

a'RZ' = -7 + ^, (461) 



a'R-^ = -5 + 3C(2) + ^ + ^L-C(2)L, (462) 

+C(2)L2, (463) 
a^R^. - -If + ^-^ + 3C(3) + 7.-6C(2).-C(3)L 

-^L' + C{2)L\ (464) 

8 Conclusions 

We have presented an exact evaluation of the master integrals with at most one 
massive propagator entering the electroweak form factor. Due to the rather simple 
threshold structure, we have succeeded in obtaining exact representations of all the 
master integrals in terms of one-dimensional harmonic polylogarithms of x. By an 
expansion of the harmonic polylogarithms for large and small values of the argument, 
we could directly obtain the large and small momentum expansions of the master 
integrals. We checked our results with previous ones present in the literature finding 
complete agreement |2Sl 1251 123 123 12H1 ■ The remaining master integrals, containing 
two or three massive propagators, have thresholds not only in s = and s = m? 
but also in s = Am? and their exact evaluation requires probably the introduction 
of new classes of special functions. 
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A Propagators 

In this appendix we list all the denominators entering the two-loop integrals com- 
puted in this paper (see section 2.2). 
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B One-loop amplitudes 

In this appendix we present the results for the one-loop diagrams which are necessary 
for the computation of the factorized two-loop amplitudes. We have recomputed 
them with the method of differential equations described in the main body of the 
paper in terms of HPL's. In some cases, we found that it was necessary to compute 
the one-loop amplitudes up to the third order in e. 
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B.l Tadpole 



^ ' i=—\ 



where: 



B.2 Bubbles 



— = -1 . 4J 
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= + (489) 



where: 



= 1, (490) 

= 2-if(0,x), (491) 

El = 4-C(2) -2i7(0,a;) + i/(0,0,x), (492) 
= 8-2(C(2) + C(3)) + (C(2)-4)if(0,x) + 2i7(0,0,a;) 

-if(0,0,0,x) , (493) 

Es = 16-4C(2)-^C'(2)-4C(3)-8if(0,a:) + 2(C(2) 

+C(3))i^(0, x) + {A- C(2))//(0, 0, a;) - 2H{0, 0, 0, x) 
+if(0,0,0,0,x) , (494) 

The above amphtude — the simplest massless one — is easily computed with a 
Feynman parameter or with the general formula in the appendix of [Zj. 



= (^)'Y,e'F, + 0{,'), (495) 
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The above amphtude is not a master integral. 
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The above amphtude does not have infrared 1/e poles because a non- vanishing 
mass in the boson line acts as an infrared cutoff on the transverse momenta: kt > m 
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C Independent two-loop amplitudes 

In this appendix we give the expressions of the reducible diagrams of figures El 
and El 



i=-2 



where: 



R 



(1) 



r'^H = --H{0,x) 



2 

5 



(1) 



R\ 



(1) 



19 

Y 

65 



- C{2) - 5H{0,x) + 2H{0,0,x) , 
-5C(2)-5C(3)-[19-2C(2)]iJ(0,x) 



(520) 
(521) 

(522) 
(523) 
(524) 

(525) 



R. 



(1) 



+ 10H{0,0,x) -AH{0,0,0,x) , 
211 11C'^(2) 

— - 19C(2) - — ^ - 25C(3) - [65 - 10C(2) 

-10C(3)]i/(0, x) + [38 - 4C(2)]i7(0, 0, x) - 20H{0, 0, 0, x) 
+8/7(0, 0, 0, 0, x) , (526) 



j=-2 



(527) 
(528) 



58 



1 

2' 
19 



-2C(2)-3i/(0,a;)-2 



1 + 



X 



H{-l,x) + H{0,0,x) 



+ 
65 



1 - 



1 



^(0,-l,x), 



(529) 
(530) 

(531) 



4C(2)-C(3)-[7-C(2)]//(0,x)-12 



1 + 



x 



H{-l,x) 



+3H{0,0,x) - 
-H{0,0,0,x) + 



1 + 



5 



X 



H{0,-l,x) + 8 



1 + 



X 



H{-l,-l,x) 



1 - 



1 



XJ 



[H{Q, 0, -1, x) - 4:H{0, -1, -1, x)] , (532) 



^ - 2C(2) - - 5C(3) - [15 - 3C(2) - 2C(3)]//(0, x) 



■[50 + 4C(2)] 
(19 + 2C(2)) 

X 



10 

X. 



H{-1, x) + [7-C(2)]i/(0, 0, x) - 17-2C(2) 



H{0,-l,x) + 48 



1 

1 + - 

a: 



i7(-l,-l,x) 



-3H{0,0,0,x) - 



1 + ^ 



[//(O, 0, -1, x) - 4/^(0, -1, -1, x) 



+ 



l + -[12i/(-l, 0,-1, x) -32if(-l, -1,-1, x)] + 



[H{0, 0, 0, -1, x) - 6H(0, -1, 0, -1, x) 

+16//(0, -1, -1, -1, x) - 4i/(0, 0, -1, -1, x)] 
+i/(0, 0,0,0,x), 



(533) 



(534) 
(535) 



1 

2' 

5 



1 

1 + - 

X. 



2 
19 

f + C(2) 



(536) 

H{-l,x), (537) 
[5//(-l, x) + H(0, -1, x) - 4i/(-l, -1, x)] 

59 



1 + 1 



-H{0,-l,x), 



(538) 



where: 



where: 



- + 5C(2)-C(3)-[19 + 2C(2)] 



1 

1 + - 

X. 



-10 



1 + 



[H{0, -1, - 2i/(-l, -1, x)] + 5H{0, -1, x) 



1 + 



[i/(0, 0, -1, x) - 4i/(0, -1, -1, x) - 6H{-1, 0, -1, x) 



+16//(-l, -1, -1, x)] + 4/7(0, -1, -1, x) - H{0, 0, -1, x) , (539) 



8 1^12 



d(4) _ 
K_2 — 


1 

2 ' 




= 


3 

2' 




d(4) _ 

— 


^ + C(2), 






y + 3C(2)- 


C(3), 




f + 7C(2) + 


9(2(2) 
5 



- 3C(3) , 



8-^13 



p(5) _ 
K_2 — 



(5) 



R\ 



(5) 



(5) 



1 

2' 
1 

2 ' 

■^-C(2), 
~-C(2) + C(3), 
i-C(2)-^ + C(3), 



(540) 

(541) 

(542) 
(543) 
(544) 
(545) 
(546) 

(547) 
(548) 

(549) 
(550) 
(551) 
(552) 
(553) 



60 



where: 



p(6) 

K_2 — 


1 

2 ' 






3 

2' 




JXq — 


^ + C(2), 






y + 3C(2)- 


C(3), 


d(6) 
^2 


y + 7C(2) + 


9C'(2) 
5 



ox2e 2 



where: 



(554) 
(555) 



(556) 
(557) 
(558) 
(559) 

3C(3) . (560) 



<r> = ,'<-./{rt.H.%}^3^^^ (561) 



(562) 



i?2 = 1, (563) 

i?5 = 3-//(0,x), (564) 

= 7-C(2)-3//(0,x) + //(0,0,x), (565) 

= 15-3C(2)-2C(3)-[7-C(2)]//(0,x) 

+3i/(0,0,x)-//(0,0,0,x), (566) 

= 31 - 7C(2) - - 2C(3) - [15 - 3C(2) - 2C(3)]if(0, x) 

+7H{Q, 0, x) - 3H{Q, 0, 0, x) + H{Q, 0, 0, 0, x) , (567) 



/.^(-^) / (568) 

^ ^ i=-2 




61 



1 

"2' 
5 

-2 + 



1 + i 

X 



H{-l,x) + -H{0,-l,x), 

X 



(570) 
(571) 



19 

-f-C(2) + 



1+i 



[5i/(-l, x) - 4:H{-1, -1, x) + 3H{Q, -1, x) 



-i/(0, -1, a;) - -[4if(0, -1, -1, x) - i/(0, 0, -1, x)] , 

X 



65 



2 -5C(2) + C(3)+[(19 + 2C(2))M+- 



iJ(-l,x) 



-20 



1 + 



X 



H{-l,-l,x)+ 10+(9 + 2C(2))- H{0,-l,x) 



X 



+ 



1 + 



X 



[3H{0,0,-l,x) + 16H{^l,-l,-l,x)-6H{-l,0-l,x) 



-12H{0, -1, -1, x)] - H{0, 0, -1, x) + 4H{0, -1, -1, x) 
+-[16H{0, -1, -1, -1, x) - 6H{0, -1, 0, -1, x) 



X 



-4:H{0, 0, -1, -1, x) + H{0, 0, 0, -1, x)] . 



(572) 




i=-2 



(573) 

(574) 



1 

2' 
3 

2' 

^ + C(2) + 



1 + i 



i/(-l,x) 



-H{0,-l,x), 



(575) 
(576) 

(577) 



-2+3C(2)-C(3) + 



1 

1 + - 

X 



[7H{-l,x) -4:H{-l,-l,x) 



+H{0, -1, x)] - 3H{0, -1, x) + 4//(0, -1, -1, x) 
-H(0,0,-l,x), 

_|,5C(2)-3C(3) + ^ 
-28 (l + -)Hi-l, -l,x) + [2 - 2C(2) + -]h{0, -1,x) 



(578) 



(33 + 2C(2))(l + -)]i/(-l,x) 



62 



+ 

+ 



1 + - 

X. 

1 

2 - - 

X. 



[16i/(-l, -1, -1, x) - 6H{-1, 0, -1, x)] 



[4i7(0, -1, -1, x) - H{0, 0, -1, x)] 

-16H{0, -1,-1, -1, x) + 6i/(0, -1, 0, -1, x) 
+4:H{0, 0, -1, -1, x) - H{Q, 0, 0, -1, x) . 



(579) 




where: 



where: 



(580) 

+ (581) 

\ ^ y ,=_9 



1 

~2 ' 

-^ + i/(0,a:), 
19 

-- + C(2) + 5HiO,x)-2HiO,0,x), 



-- + 5[C(2) + C(3)] + [19 - 2C(2)J H{0, x) - 10H{0, 0, x) 



(582) 
(583) 
(584) 

(585) 



65 

y 

+4//(0, 0, 0, x) , 

211 ll(r^f2) r 

-— + 19C(2) + — ^ + 25C(3) + [65 - 10(C(2) 

+C(3)) H{0,x)- 38-4C(2) H{0,0,x) + 20H{0,0,0,x) 
-8H{0,0,0,0,x). (586) 




^ ^ i=-3 



(587) 
(588) 





1 

a; 






(589) 






-H{0, 


^)], 


(590) 


R^l? = 




-C(2) 


-i/(0,a;) + 1^(0,0, x)], 


(591) 




-> 


-C(2) 


- 2C(3) - (1 - C(2))//(0, x) + //(O, 0, x) 





63 



-i/(0,0,0,x)], 



R'i'^ = -- 1-C(2)-2C(3) 



1 r 



9(2(2) 



(l-C(2)-2C(3))//(0,x) 



XL 10 

+ (1 - C{'^))H{0, 0, x) - H{0, 0, 0, x) + H{0, 0, 0, 0, x) 



where: 



(12) 



6C(3) 



x 



^C'(2) + 12C(3)-12C(3)if(0,x) 



1 



i=0 



(592) 

(593) 

(594) 
(595) 

(596) 

(597) 

(598) 
(599) 



where: 

aR 



(13) 



aRl 



(13) 



= -[a2)H{-l,x) 

X 

+2H{-1, 0, -1, x) + i/(-l, 0, 0, x) + 2H{0, -1, -1, x) 
-H{0,-l,0,x)], (600) 

= i[2(C(2) -2C(3))//(-l, x)-2C{2)H{0-l, x)+C{2)H{-l, 0, x)) 

-2C{2)H{-1, -1, x) - 4//(-l, 0, -1, x) + 2//(-l, 0, 0, x) 
+4:H{0, -1, -1, x) - 2H{0, -1, 0, x) + 4:H{-1, -1, 0, -1, x) 
-2//(-l,-l, 0, 0, x)+4:H{-l, 0,-1,-1, x) +2H{-1, 0,-1, 0, x) 
-2H{-1, 0, 0, -1, x) -3//(-l, 0, 0, 0, x)-12H{0, -1, -1, -1, x) 
+2H{0, -1, -1, 0, x) + 4:H{0, -1, 0, -1, x) + 2H{0, -1, 0, 0, x) 
+2H(0, 0, -1, -1, x) - //(O, 0, -1, 0, x)] , (601) 

= (^)"E^'^1"^ + ^(^')' (603) 



64 



(14) 



alt. 



(14) 



aR\ 



(14) 



1 

2 ' 
1 

2 ' 

^ + C(2), 

i + C(2) - C(3) , 




(15) 



02) 
2x ' 

^[C(3)-2C(2)i/(0,x)], 



1 



X 



e{2) 

10 



+ C(3)i/(0,a;)-2C(2)i/(0,0,a;) 




^ ^ i=-3 



[l-^(0,x)] 



= -[2-C{2)-H{0,x)-(i + -)h{-1,x) + H{0,0,x) 

+H{0,-l,x)], 
-- -U- C(2) - 2C(3) - (1 - C(2))i/(0, x)-3(l + -)h{- 

XL \ X/ 

+H{0, 0,x) + a(i + H{-1, -1, x) - H{0, -1, x) 



65 



-H{0,0,0,x) + H{0,0,-l,x) -4:H{0,-l,-l,x) , (618) 
1 Tg + ^(2) - - 2C(3) - (1 - C(2) - 2C(3))i/(0, x) 



X 



-(7 + 2C(2)) (l + i) H{-1, x) + (l- C(2))//(0, 0, x) - (5 
-2C(2))//(0, -1, x) + 12 (l + i/(-l, -1, x) - H{0, 0, 0, x) 

-//(O, 0, -1, x) + 4//(0, -1, -1, x) + (l + [6//(-l, 0, -1, x) 

-16H{-1, -1, -1, x)] + iJ(0, 0, 0, 0, x) + i/(0, 0, 0, -1, x) 
+16H{0, -1, -1, -1, x) - 6H{Q, -1, 0, -1, x) 

-4i/(0,0,-l,-l,a;)l , (619) 



,H'-n^Ji,W>=,}^^^^^ (620, 

^ ^ i=-l 




--i7(l,0,x), (622) 
--[2C(2)//(1, x) + 2H{1, 0, x) + //(0, 1, 0, x) 

X 

-2//(l,0,0,x) + 2//(l,l,0,x)], (623) 

-i[4(C(2) + C(3))i/(l,a;) + 2C(2)i7(0,l,x) + 2(2-C(2))i/(l,0,a;) 

X 

+4C(2)i7(l, 1, x) + 2H{0, 1, 0, x) - AH{1, 0, 0, x) 

+4:H{1, l,0,x) + H{0, 0, 1, 0, x) - 2H{0, 1, 0, 0, x) 

+2H{0, 1, 1, 0, x) + AH{1, 0, 0, 0, x) + 2if(l, 0, 1, 0, x) 

-4H{1, 1, 0, 0, x) + AH{1, 1, 1, 0, x)] , (624) 




(626) 



(627) 



-1-[2H{-1, x) - AH{-1, -1, x) + H{0, -1, x)] , (628) 

- — [2(2 + ((2))H(~-1, x) + 2H(0. ~l.x) - 8H(-1. -l,x) 
2x 

+ 16H{-l,-l,-l,x) -6H{-l,0,-l,x) - 4:H{0,-l,-l,x) 
+H{0,0,-l,x)], (629) 



(^j^E/'-sr+fi''). (631) 




i|l-/f(0,i)], (633) 

i [2 - C(2) - 2H{Q, x) + 2i/ (0, 0, x)] , (634) 
i[4 - 2C(2) - 5C(3) - 2(2 - C(2))i/(0, x) + 4/7(0, 0, x) 

X 

-4i7(0,0,0,x)], (635) 
- [S - 4C(2) - - 10C(3) - 2(4 - 2C(2) - 5C(3))i/(0, x) 

XL 5 

+4(2-C(2))//(0, 0, x) -8//(0, 0, 0, x)+8H{0, 0, 0, 0, x)], (636) 



= (^)"E^^^f°^+^(^')' (638) 




--//(l,x), (639) 

--[2H{-l,x) -2H{-l,-l,x) - H{-l,0,x) 

+H{Q,-l,x)], (640) 

i[-4i/(-l, x) + AH(-1, -1, x) + 2H(-1, 0, x) - 2//(0, -1, x) 

a; 

-4i/(-l, -1, -1, x) - 2i/(-l, -1, 0, x) + 2H{-1, 0, -1, x) 
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aR\ 



(20) 



-2H{-1, 0, 0, x) + 2H{0, -1, -1, x) + H{0, -1, 0, x) 
-H{0,0,-l,x)], (641) 

^[(8 - 6C(3))//(-l, x) + 4i/(0, -1, x) + 2C(2)//(-l, 0, x) 

-4:H{-1, 0, x) - 8i/(-l, -1, x) + 4i/(-l, 0, 0, x) 
-AH{-1, 0, -1, x) - AH(0, -1, -1, x) + 2H{0, 0, -1, x) 
-2H{0, -1, 0, x) + 4//(-l, -1, 0, x) + 8H{-1, -1, -1, x) 
-8H{-1, -1, -1, -1, x) - 4:H{-1, -1, -1, 0, x) 
+4//(-l, -1, 0, -1, x) - 4//(-l, -1, 0, 0, x) 
+AH{-1, 0, -1, -1, x) + 2//(-l, 0, -1, 0, x) 
-2H{-1, 0, 0, -1, x) - 4//(-l, 0, 0, 0, x) 
+AH{0, -1, -1, -1, x) + 2H{0, -1, -1, 0, x) 
-2H{0,~l,0,-l,x) + 2H{0, -1,0,0, x) 
-2H{0, 0, -1, -1, x) - H{0, 0, -1, 0, x) 

+H{0,0,0,-l,x)], (642) 



where: 




'--2 



(21) 



2(4-D) /l^^i^f^ll^i^i^f)) ^ (643) 

V / i— _9 



i/(0,-l,x), 

[if(0, 0, -1, a;) - AH{0, -1, -1, x)] 



2x 
1 

2a; 

^[2C(2)i/(0, -1, x) + i/(0, 0, 0, -1, x) - 4//(0, 0, -1, -1, x) 
-6i/(0, -1, 0, -1, x) + 16i/(0, -1, -1, -1, x)] , 



(645) 
(646) 

(647) 



where: 




aR 



(22) 



(648) 

(649) 



4x ' 



(650) 
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where: 



where: 







(651) 




-^[C(2)- 2/7(0, 0,0;)], 


(652) 




-^[5C(3) - 2C(2)/7(0, x) + 4/7(0, 0, 0, x)] , 

_ 1 rll^ _ 5^^3^^^Q^ ^ 2C(2)//(0, 0, x) 
xl 10 


(653) 



-4//(0,0,0,0,x)], 




-//(0,0,0,x) + //(0,0,0,0,x)], 




4X>5l>iO 



a/? 



(24) 



(24) 



a/? 



^[l-//(0,a;)], 

-[2-C(2)-2//(0,x) + 2//(0,0,x)], 



(654) 

(655) 
(656) 



aR^S = 
aR^lf = 


1 


-H{0,x)], 


(657) 
(658) 


aR^^f = 




H{0,x) + H{0,Q,x)], 


(659) 


aR^^ = 
a/?f^) = 


>- 

4 r 

- 1 - 
X - 


C(3) - H{0, x) + H{0, 0, x) - H{0, 0, 0, x)] , 
- C(3) - - (1 - C(3))//(0, x) + //(O, 0, x) 


(660) 



(661) 

(662) 
(663) 



(664) 
(665) 
(666) 
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-[2 - C(2) - C(3) - (2 - C(2))i/(0, x) + 2H{0, 0, x) 

X 

-2H{0,0,0,x)], (667) 
4 r C^(2) 

- 4 - 2C(2) - ^ - 2C(3) - 2(2 - C(2) - C(3))i/(0, x) 
+2(2-C(2))i/(0,0, x) - 4//(0,0,0, x) + 4i/(0,0,0,0, x)] , (668) 



(669) 
(670) 



-ii/(-l,x), (671) 

X 

--[2H{-1, x) + H{-1, 0, x) - 2H{-1, -1, x)] , (672) 

-^[(4 - a2))H{-l,x) - 2H{-l,0,x) - 4H{-1, -l,x) 
+4if(-l,-l,-l,x) +2H{-l,-l,0,x) + H{-l,0,0,x)], (673) 
-i{2(4 - C(2) - CmH{-l, x)-{4- amH{-l, 0, x) 

-2(4 - C{2))H{-1, -1, x)] + 2H{-1, 0, 0, x) 

+4:H{-l,-l,0,x) + 8H{-l,-l,-l,x) - H{-l,0,0,0,x) 
-2H{^1, -L 0, 0, x) - AH{-1, -1, -1, 0, x) 
-8H{-l,-l,-l,-l,x)}, (674) 



= .^^-^^{^^^m^^^.}^;^^ (675) 
- (^)"E^^^r^+^(^')' (676) 

^ ' i=— 3 




(677) 

^[1-H{0,x)], (678) 
-[1 - C(2) - H{0, x) + H{0, 0, x)] , (679) 
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= l[l-C(2)-2C(3)-(l-C(2))i/(0,a:) + //(0,0,x) 



X 



-H{0,0,0,x)], (680) 



= ^ i_^(2)-2C(3)-:^-(l-C(2)-2C(3))/f(0,x 



1 



X 



9(2(2) 



10 



+ (1 - C(2))i^(0, 0, x) - H{0, 0, 0, x) + H{0, 0, 0, 0, x)] . (681) 
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